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PAUL BAIRD 

Abstract. We study a class of complex polynomial equations on a fi- 
nite graph with a view to understanding how holistic phenomena emerge 
from combinatorial structure. Particular solutions arise from orthogo- 
nal projections of regular polytopes, invariant frameworks and cyclic 
sequences. A set of discrete parameters for which there exist non-trivial 
solutions leads to the construction of a polynomial invariant and the 
notion of a geometric spectrum. Geometry then emerges, notably di- 
mension, distance and curvature, from purely combinatorial properties 
of the graph. 
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1. Introduction 

Given a graph T = (y,E), with vertices V and edges E together with a 
function ip : V ^ C, the equations we wish to study ah have the form 

^ ' yr^x yr^x 

at each vertex x, where y ~ x means y is adjacent to x, n(x) is the degree 
of r at X (the number of vertices adjacent to x) and where 7 : y — t- R is 
a rea/-valued function, which in many situations we will suppose constant; 
then we call the solution regular. A special case is when 7 = and we have 
the equation: 

(2) = 5](v9(x)-(^(y))^ 

y~x 

for all X € V, solutions of which we refer to as holomorphic functions. Note 
that the equations are invariant by the replacement of 99 by = Xip + /i 
where A and /i are complex numbers, as well as with respect to complex 
conjugation ip ^Tp. 

In what follows, it is convenient to write Aip{x) = Ylyr^x {'P(.y)~^{^)) 
(the Laplacian) and {d(p)'^{x) = ■:^'}2yr^x (fiv) ~ V^i^))'^ (t^e symmetic 
square of the derivative), whereby equation ([1]) becomes: 

7(x)A93(x)^ = (d(^)^(x) . 

All our graphs will be simple, that is we do not allow loops or multiple edges, 
although much of what we discuss can be generalized to non-simple graphs. 

For a given graph T = {V,E), we are particularly interested in the set S 
of constant values of 7 that can occur for which ([1]) has non-trivial solutions: 

S := {7 S M : 3 non — constant (f : V ^ C such that jAip'^ = (d(^)^} . 

We will call this set the geometric spectrum of T, which we regard as an 
intrinsic object associated to the graph, somewhat akin to the spectrum of 
the Laplacian. We will see that this "spectrum" reflects geometric proper- 
ties, as well as combinatorial properties of T. Note that on a finite graph, 
by the maximum principle, ip non-constant implies Aip ^ 0, so that T, is 
always well-defined. 

Particle physics provides some motivation to study these equations. As 
a lattice model, equation ([T]) expresses the influence exerted on a particular 
vertex by each of its neighbours. In the case of a graph of degree three, in 
[3] it is shown that equation ^ provides a discrete analogue of the equation 
for a shear-free ray congruence on space-time. In the smooth setting such 
congruences can be used to generate solutions to the zero rest-mass field 
equations. In Section [9l we refer to a connected graph F as a particle and 
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an equivalence class of solutions to ([T]) on F as a state of that particle. States 
which have 7 constant will be referred to as isostates. These are significant, 
since empirical evidence suggests that isostates occur as extremal states of 
a natural energy functional derived in Section [UJ We take the view that (i) 
at a fundamental level, the description of the world should be combinatorial 
and (ii) the only physical quantities that have meaning are relative. 

The idea that combinatorial structures should be at the basis of the de- 
scription of matter, was put forward by R. Penrose in a well-known paper of 
1971 [23]. In this article, he introduced the idea of a spin network, which is 
a graph of degree three whose edges are labeled by integers which represent 
twice the angular momentum. The notion has been generalized in recent 
years by C. Rovelli and L. Smolin in their development of loop quantum 
gravity, see [26]. Another key idea of Rovelli, is the relational interpretation 
of quantum mechanics, which puts forward the thesis that the only quanti- 
ties that have meaning are relational. Thus no particle exists in an absolute 
state; the only meaningful way to describe a particle is how it correlates 
with other particles, or systems |25j . 

This philosophy is consistent with the idea that an isolated system be 
described by a graph coupled to a complex field. Firstly, a graph is a com- 
binatorial structure that expresses a binary relation between its vertices. 
We do not need to imagine the graph as lying in some ambient space; the 
only relation that matters is whether one vertex is joined to another or not. 
Of course, one may envisage more complicated combinatorial relations in 
nature, so perhaps edges should be labeled by an integer, or by a group rep- 
resentation, as in loop quantum gravity. Secondly, the introduction of the 
field (/? on the graph, represents an additional structure from which geometry 
emerges. The equations which govern this field should enjoy the invariance 
(p I—)- Xip + fj,, so that only relative values of (p are significant. 

The condition that 7 be real in equation ([T|) is justified by Theorem 16.11 
This expresses the fact that under reasonable hypotheses, there is a unique 
realization of the solution about a vertex as a regular star in a Euclidean 
space; this enables us to develop geometry, in particular curvature and dis- 
tance. The constancy of 7 then implicates global properties of the graph: it 
is a uniformizing parameter. It is the case that the framework of a regular 
polytope, when projected orthogonally into the complex plane, determines 
a solution to ([T|) with 7 constant (Theorem 14. 4p . 

The factor l/n{x) on the left-hand side of ([T|) is explained by equation 
(|12|) of Corollary 14. 2t this means that 7 depends only on local geometric 
invariants and not on the degree. This choice also provides a consistent 
boundary, once more independent of the degree, between real solutions and 
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complex solutions (see Lemma l3.6p : there are also heuristic arguments that 
arise from the functional analytic perspective given in Appendix ICl 

A notion of holomorphic function somewhat similar to ours has been 
introduced by S. Barre However, in addition to Barre requires that 
(p be harmonic. This is quite restrictive and in particular, by the maximum 
principle, with this definition no non-constant holomorphic function could 
exist on a finite graph. 

We begin our paper with some elementary remarks concerning graph 
colourings. The next three sections then provide examples of solutions to 
equation ([1]) which will serve as a reference for subsequent developments. 
We give a complete characterization on a cyclic graph; for 7 constant, real 
solutions can be described in terms of polynomial equations with integer co- 
efficients (Theorem 13. ip . Complex solutions correspond to polygonal chains 
(or planar bar frameworks) whose edges all have the same length. In Ap- 
pendix [B1 we characterize the regular polygons or star polygons as those 
solutions having 7 constant. Indeed the energy functional that we derive in 
Section [6] (see below) defines a gradient fiow on the configuration space of 
polygonal chains which we conjecture determines an evolution towards the 
most compact regular configuration. 

Following work of Eastwood and Penrose [T^, in Section[3]we discuss how 
solutions arise from orthogonal projections of regular polytopes. An impor- 
tant step is Corollary 14.21 which shows that a certain class of star graphs in 
M.^ have invariant properties. In particular, their projections to C satisfy 
([1]) at the central vertex with 7 independent of the position of the star. By 
extending the class of invariant star graphs in Section [5l we describe a more 
general family of invariant frameworks which define solutions to ([1]). For 
example, we show that it is possible to place a complete graph on any num- 
ber of vertices in in an invariant way (Corollarv 15. 4p . By this we mean 
that the projections to C of the vertices satisfy ([T]) with 7 constant, indepen- 
dently of any rigid motion or dilation of the framework. More sophisticated 
invariant structures are described in Appendix [Al When we change our per- 
spective and begin with a combinatorial structure together with a spectral 
value 7, these examples take on added significance as "geometric objects" . 

In Section [6] we consider the converse question of how to lift a solution 
to (dJ at each vertex to an embedded star. This will enable us to associate 
both distance and curvature to a graph P coupled to a field tp, by considering 
at each vertex the best-fit polytope with regular vertex figure. We call this 
the lifting problem; its solution, expressed by Theorem 16.11 is at the heart 
of our study. It transpires that when we lift to dimension three and only to 
this dimension, the solution is the unique minimum of a natural functional 
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determined by the field i-p. This energy functional plays an important role 
in the formulation of our elementary universe in Section [9] (Definition I9.3p . 

Curvature is defined in Section [7] by analogy with the well-known notion 
for polytopes embedded in a Euclidean space, in terms of angular deficiency. 
A different, edge- curvature^ measures the angular tilt between the axes of 
adjacent vertex figures. This leads to analogues in a discrete setting of 
sectional and Ricci curvatures. 

An important construction is that of a polynomial invariant associated to 
a finite connected graph F. This is defined to be the least degree univariate 
polynomial p-^ in the spectral parameter 7 lying in the ideal generated by 
the equations ([TJ with 7 now constant and complex, after we have taken into 
account the normalization 99 1— )• + [i. The geometric spectrum occurs as 
real roots of this polynomial, however, there may be other real roots that 
are not in the spectrum. In Section [HI we use Grobner bases to find "p^ for 
a number of examples. 

Our ultimate objective is to describe an elementary universe populated 
entirely by graphs, from which geometry and dynamics emerge. For an indi- 
vidual graph, its geometry is implicit in the geometric spectrum. However, 
the nature of this geometry (which spectral value applies) should only be- 
come manifest upon correlation with another graph: it is a relative concept. 
Dynamics corresponds to change which occurs when two graphs combine, a 
graph mutates, or vertices and edges are created or annihilated. With an 
appropriate definition of thermal time, the universe is then endowed with 
local geometry and time. 

2. Elementary observations on graph colourings 

If a vertex x has degree one, that is, it has precisely one neighbour y, 
then if is any function with (/?(y) 7^ "pix)-, equation ([1]) is satisfied at x, 
with 7(x) = 1. More generally, if we can colour the vertices of the graph 
with two colours in such a way that every vertex is adjacent to at most one 
vertex of a different colour, then once more ([T]) is satisfied with 7(x) = n(x) 
for all X . To see this, we just associate the values and 1 to the two 
colours. This is easily generalised to the following. 

Lemma 2.1. Let F = {V^E) he a graph whose vertices are assigned two 
colours, and 1 say, and let ip{x) = colour of vertex x, for each x £ V. 
Suppose F is coloured in such a way that for some function k{x), either each 
vertex x is adjacent to precisely k[x) of a different colour, or each vertex 
is adjacent to no vertex of a different colour. Then ip satisfies ([T]) with 
7(2;) = n(x)/k{x) for all x £V. 
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So solutions to ([T]) arise from particular 2-colourings of the graph. But 
what about 3-colourings, or m-colourings? 

Consider the three colours {0, 1, | + ^i}, corresponding to the positions 
of the vertices of an equilateral triangle in the plane. Now, whatever vertex 
X we take, if it is adjacent to precisely one of each of the other two values, 
equation ^ is satisfied with 7 = n(x)/3. Similarly, a colouring using four 
colours with this property (with 7 replaced by n{x)/4) is given by the four 
complex numbers {0, 1, i, 1 + i}. These points correspond to the projections 
of the vertices of a regular tetrahedron in onto the complex plane. 

A result of Eastwood and Penrose [T5], affirms that the coordinates of 
any orthogonal projection of the vertices of a regular simplex in into C, 
provide + 1 colours that can be used to generate solutions to (P). Specifi- 
cally, let {vi,V2, . . . , WAf+i} be the vertices of a regular A^-simplex in and 
let P : — )• C be any orthogonal projection. Let S := {zi, Z2, ■ ■ ■ , zn+i}, 
where Zk = P{vk) (A; = 1, . . . , iV + 1). 

Lemma 2.2. Let T = [V^ E) he a graph whose vertices are assigned A'' + 1 
colours taken from the set S. Let (p{x) = colour of vertex x. Suppose T is 
coloured in such a way that either each vertex is adjacent to precisely one 
of each of the other N colours, or each vertex is adjacent to no vertex of 
a different colour. Then ip satisfies ([T]) with 7(x) = n{x)/{N + 1) for all 
xeV. 

Example 2.3. A bipartite graph Km,n is a graph with m + n vertices which 
can be separated into two sets: and with the 

property that for each j = 1, . . . , m, we have Xj ~ yr for all r = 1, . . . , n; 
Xj 9^ Xfc for all j ^ k; yr yg for all r ^ s. The bipartite graph Kn^n 
admits solutions deriving from both of the above lemmas. In the first case, 
we can colour the vertices xi, . . . ,Xn with and the vertices yi, . . . ,yn with 
1. Then the corresponding function ip solves ([T]) with 7 = 1. In the second 
case, we let both xi, . . . ,Xn and yi, . . . , y„ to be the orthogonal projections 
zi, . . . , Zn of the vertices of a regular (n — l)-simplex in M"""*^, in any order. 
Once more ([1]) is satisfied with 7 = 1. 



yi y2 ys 




Xi X2 X3 



We can explicitly find all solutions to ([T|) on the bipartite graph 
illustrated above when 7 = 1. If we first normalize so that xi = and 
yi = 1, then the solutions are given by a;2 = X3 = 0, 7/3 = A an arbitrary real 
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parameter, y2 = ^^^-^ + A ^ ^^2^ ) • ^ varies from 1 to 0, the solution 

interpolates between 2/1 = 2/2=2/3 = ! and 2/1 = Ij 2/2 = ^^5^, 2/3 = 0, which 
are the vertices of an equilateral triangle. 

If we now normalize so that xi = and X2 = 1, then necessarily X3 = 
^^2^ ' then have a 2-parameter family of solutions given by arbitrarily 
prescribing 2/2 = A, 2/3 = and then 2/1 = A ^ ^^2^ ) + ( ^^2^ ) • the 
case when ^ = and A = 1, we once more have the vertices of an equilateral 
triangle, corresponding to the hypotheses of Lemma 12.21 It turns out that 
for ivr33, 7 = 1 is the only possible constant value of 7 for which ([1]) has 
non-constant solutions. This will be justified in Section [8l 



3. Regular cyclic sequences 

We shall call a solution to ([T|) on a cyclic graph with 7 constant, a regular 
cyclic sequence. Real regular cyclic sequences have a particularly simple con- 
struction from polynomial equations with integer coefficients, as we explain 
below. First we deal with the trivial case when 7 = 2. 

Let (xo, xi, X2, . . . , xat-i, xat = xq) be a regular cyclic sequence which 
solves ([1]) with 7 = 2. Consider a particular segment of three successive 
terms Xfc_i, Xfc, Xjt+i. On applying ([1]) at the vertex Xfc, we obtain the equa- 
tion: 

(xfc+i + Xk-i - 2xkf = (xfe-i - Xkf + (xfc+i - Xkf 
<^ {xk - Xk-i){xk - Xk+i) = , 

so that necessarily, Xk is equal to one of its neighbours. Conversely, as 
described in the previous section, if every vertex is adjacent to at most one 
of a different value, then the sequence solves ([1]) with 7 = 2. Thus any cyclic 
graph of order > 4 admits such sequences: the cyclic graph is coloured by 
numbers (real or complex) in such a way that at each vertex, at least one of 
its neighbours carries the same colour. This amounts to colouring the graph 
so that connected segments of order at least two have the same colour. We 
now consider the general case. In what follows, we refer to normalization as 
the freedom ^ Xif + n (A,// G C). By a real regular cyclic sequence, we 
mean one in which every term is real under some normalization. 

Theorem 3.1. (Construction of real regular cyclic sequences).' Take any 
polynomial a„x"-|-a,i_ix"~^+- • •+aix+ao with integer coefficients all strictly 
positive. Multiply the polynomial by x + 1 to obtain the new polynomial 

p{x) := 6„+ix"+^ + 6„x" + --- + bix + bo 

= Onx""^^ + (a„ + a„_i)x" H \- {ai + ao)x + oq • 
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Let X = y be any real root of p(x). Then a cyclic sequence 

{xo,Xi,X2, ■ ■ ■ ,XN-1,XN = Xo) 

of order N = 2 is constructed by arbitrarily prescribing xq and then 

requiring increments yi = Xi — Xi-i of successive terms to be taken from the 
set {1, y, y^, . . . , y""^^} in such a way that each increment y^ occurs precisely 
bk times and any two adjacent increments have powers that differ by precisely 
one. This is always possible and up to these constraints, the ordering is 
arbitrary. The constant j in is given by j = 2(1 + y^)/(l — y)^. 

Conversely, up to a multiple, addition of a constant and cyclic permu- 
tations, any real regular cyclic sequence with 7 7^ 2 or 1 arises this way. 
The real regular cyclic sequences with 7 = 2 are characterized as those 
made up of connected segments of order > 2 on which the sequence is con- 
stant; those with 7 = 1 oscillate and up to normalization are equivalent to 
(0,1,0,1,... ,0,1). 

We refer to the increment y in the above theorem, as the fundamental 
increment associated to the real cyclic sequence. 

Remark 3.2. Since any root y must be strictly negative and adjacent pow- 
ers differ by one, it follows that a real regular cyclic sequence must oscillate. 
The length of the sequence is given by Y2k^k = 2^^afc, so that a non- 
trivial sequence can only occur on a cyclic graph of even order (which is also 
a consequence of oscillation) . 

Example 3.3. If we take for our starting polynomial x -\- 2, then multipli- 
cation by X + 1 gives the polynomial + 3x + 2 with real root x = —2. 
We can now arrange the powers of this root with appropriate multipicity to 
give the sequence of increments (1, y, 1, y, y'^,y) = (1, —2, 1, —2, 4, —2). We 
construct a real regular cyclic sequence of order 6 by first setting xq = 
and then proceeding so that xi — xq = 1,X2 — x\ = —2 and so on. We 
thereby obtain the sequence (0, 1, —1, 0, —2, 2) on a cyclic graph of order 6. 
Since the sequence is only defined up to multiple, addition of a constant and 
cyclic permutations, we can normalize the sequence in such a way that the 
minimum value is and that this occurs for the first term: (0, 4, 2, 3, 1, 2). 



1 3 
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Example 3.4. Irrational sequences arise from irrational roots. For example, 
let us start with the polynomial + 4x + 1, with root x = —2 + ^/?,. On 
multiplying by x + 1 we obtain the polynomial x'^ + hx'^ + 5x + 1 . A suitable 
sequence of increments is given by {l,y,y^,y,y'^,y,y'^,y,y^,y'^,y^,y) with 
y = —2 + \/3. On calculating, we can now construct a real regular cyclic 
sequence of order 12; explicitly, it is given by (0, 1, -1 + \/3, 6 - 3\/3, 4 - 
2V3, 11 - eVs, 9 - bVs, 16 - 9V3, -10 + 6^/3, -3 + 2^/3, -5 + SVS, 2-^/3). 
Since the absolute value of the root is < 1, all terms of this oscillating 
sequence lie in the interval [0, 1]. The value of the constant 7 in ([TJ is given 
by 7 = 4/3. 

In order to prove Theorem 13. II we first of all establish a recurrence relation 
that determines a subsequent term of the sequence in terms of three previous 
terms. 

Lemma 3.5. Let (xo,xi, . . . ,xn-i,xn = 2:0) 0, non-constant regular 
cyclic sequence satisfying ^ with 7/2, then the increments yt = Xk+i — Xk 
satisfy the recurrence relation: 

either yk-i^/yk-2 
or yk-2 ■ 

Conversely, any sequence of increments satisfying these relations determines 
a real regular cyclic sequence. 

Proof. Let {xo,xi, . . . ,X]y-i,XN = xq) be a non-constant regular cyclic se- 
quence satisfying ([T]) with 7 7^ 2. Consider a particular segment of the 
sequence consisting of four consecutive vertices: {xk-2,Xk-i,Xk,Xk+i)- On 
normalising, we can suppose this segment equivalent to {k,0,x,y), for some 
real numbers k, x, y. Note that k ^ 0, for otherwise x would have to be 
zero since if not, we would have 7 = 2 at vertex k — 1. But then proceed- 
ing along the cycle, we would eventually encounter a non-zero value at a 
vertex, which would then imply 7 = 2 at that vertex; a contradiction. On 
evaluating equation ([1]) at the vertex Xk-i, we obtain: 

^ 2{k^ + x^) 

Now evaluate ([T]) at vertex x^: 

7(y-2x)2 = 2((y-x)2 + x2). 

On eliminating 7, we obtain the quadratic equation in y: 

ky"^ + {k — x)'^y — x{k — x)^ = . 

This gives the two possible values x{k — x)/k and x — k for y. To recover 
the general case, we set k = kj._2 — Xk~i, x = Xk — Xfe_i, y = x^+i — Xfc_i. 
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This gives the two values: 

Xk-l{Xk - Xk-l) + Xk{Xk-2 - Xk) 



Xk+1 = \ Xk-2 - Xk-l 

Xk + Xk-l - Xk-2 

On subtracting Xk from both sides, we obtain the recurrence relation for the 
increments in the statement of the lemma. The converse can be shown by 
direct computation. □ 

Proof of Theorem 13.11 Let {xq^xi, . . . ,xn-i^xn = xq) be a non-constant 
real regular cyclic sequence satisfying ([T|) with 7 7^ 2. From Lemma 13.51 
the sequence of increments (y^ = x^+i — Xk) must alternate in sign, for 
otherwise, if we have two consecutive increments of the same sign, then 
all subsequent increments would have the same sign and the sequence {xk) 
would be monotone increasing or decreasing, which is impossible. 

First normalise so that xi — xq = 1, X2 — xi = x < and consider the 
sequence of possible increments: 



l,x, < 



x^, 



1, 



x^, 



X, 

1/x, 

X, 



Suppose first that x = —1, then we obtain the sequence of increments 
(1,-1,1,-1,. ..,1,-1) which corresponds to the real regular cyclic sequence 
(0, 1,0,1,... ,0, 1) with 7=1. Furthermore, any non-constant regular cyclic 
sequence with 7 = 1 must have this form, since if we take a segment 
Xk-l ~ Xfc ~ Xfc+i, then 

{xk-Xk-i+Xk-Xk+if = 2{xk-Xk-if+2{xk-Xk+if <^ (xfc+i-Xfc-i)^ = 0, 

so that Xk+i = Xk-l- Henceforth, suppose that x 7^ —1. In particular, since 
X is real, we cannot have x*^ = 1 for any power /c 7^ 0. 

Note that every term in the sequence of increments must be a power of x. 
Now multiply through by the highest negative power of x to obtain 1 = x'^ 
in some position, with all other powers of x greater than or equal to zero: 

(x^o,x^\...,l,...), 

where each > 0. Now cyclically permute the sequence to obtain 1 in the 
first entry: 

where k is necessarily an odd positive integer. The recurrence relation then 
implies that all terms must have the form [x^'Y for some integer r > 0, so 
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we now set y = x , to obtain the sequence: 




On applying the recurrence relation once more, we see that an occurrence 
of must be followed by either y^~^^ or y^~^. But X^^Zg = Z]fc=o(^fc+i — 
Xk) = 0, which implies that y satisfies a polynomial equation of the form: 

p{x) := ar+ix""^^ + (a^+i + + 2x''~^ + 2x''~'^ H 

r 

(4) • • • + + (ao + l)x + ao + Y^ I3s{x' + x*"^) 

3 = 2 

= {x + l)(ar.+ix^ + x^~^ + x^~2 + ■ ■ ■ + X + ao) 

r 
s=2 

where ar+i,ao are strictly positive integers and fig are integers that are 
> 0. Thus for each of the a^+i occurrences of the maximum power y^"*"^, we 
must have at least one more occurrence of y*". Similarly for each of the ao 
occurrences of the minimum power y^. Since by assumption y 7^ 1 so that 
y^ ^ 1 for any k ^ 1, all intermediate powers must occur at least twice. 
However, we may have further oscillations between powers of and y^~^ 
for s = 2, . . . , r, which are given by the coefficients fSg- Then the polynomial 
p{x) has the form of the statement of the theorem. 

Conversely, given a polynomial anX^ + a„_ix"'~^ + • • • + aix + ao, as in 
the statement of the theorem, it can be written uniquely in the form: 

r 

{ar+ix'' + x""^ + x"""^ + • • • + X + ao) + ^ Psx'~'^ . 

s=2 

The expression for 7 is deduced from In that expression x and —k 
are successive increments, so we must have x/k = —y or —y~^. But the 
expression for 7 is invariant under y 1— )• y~^. □ 

We now consider complex regular cyclic sequences, that is solutions to ([1]) 
with 7 constant which in any normalization have at least one non-real value. 
We first prove a lemma in a more general context, which we will require in 
later sections. 

Lemma 3.6. Suppose the equation: 

(n \ 2 n 
1=1 / £=1 

is satisfied for 7 real and for zi real and not all zero. Then 7 > 1. In 
particular, if ip solves ^ at a vertex x of degree n, if n>2 and 7 < 1, then 
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in any normalization which has ip{x) = 0, at least one of ip{y) (y ^ x) must 
he complex. 

Proof. The Cauchy-Schwarz inequality shows that for any set {ai, . . . ,a„} 
of complex numbers, one has the inequality 

n / " \ / " \ 

(5) n ^ aiat > I ^ a J I ^ | 

with equality if and only if ai = a2 = • • • = 0^. Then for real and not all 
zero satisfying ([T|), we have: 

2 

2 



□ 



Corollary 3.7. Let (xq, xi, X2, • • • , XAr_i, xat = xq) he a regular cyclic se- 
quence. Then either there is some normalization under which every term is 
real, or whatever normalization is taken, every consecutive triple {xk-i,Xk, x^+i) 
contains at least one complex term. The real regular cyclic sequences are 
characterized by the property 7 > 1. 

Let us explore in more detail the import of the above corollary. Let 
(xo, xi, X2, . . . , XAr_i, xtv = xo) be a regular complex cyclic sequence. Then 
for each term x/j, whatever the normalization, not all of (xfc_i, x^, x^+i) are 
real. We can normalize so that x^ = and its two neighbours take values 
as indicated in the diagram: 



where we suppose r,s > and 9 G (— 7r,7r) \ {0}. On substituting into ([T]) 
we obtain 

7 = 2 + 



^2g-ie _|_ g2gie _ 2rs ' 
This is real if and only if 

2 cos 9 

(6) r = s m which case 7 = , 

cos 9 — 1 

which confirms the inequality 7 < 1. The case when 9 = ib7r/2 corresponds 
to 7 = 0, that is holomorphicity at the vertex x. We underline the fact that 
the two cases, real triple or complex triple, are mutually exclusive, their 
nature determined by the value of 7. 
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Thus the terms of any complex regular cyclic sequence determine a closed 
walk in the plane, such that at each step the walk progresses along an edge 
of fixed length in such a way that the angle between successive edges is ±9, 
for some fixed 6; the latter property being a consequence of the constancy of 
cos^ in ([6]). In particular, any regular polygon satisfies these criteria, with 
9 the (constant) exterior angle. 

In principle it is possible to find all solutions: if 9 is the exterior angle and 
we perform k rotations through +6 and i through —9, then {k — i)9 must be 
a multiple of 27r with k + i = N , the order of the graph. Thus for each A^, 
up to normalization and conjugation, there can only be a finite number of 
possibilities. In particular, taking into account the real cyclic sequences, the 
geometric spectrum of a cyclic graph is finite. The regular cyclic sequences 
on the cyclic graphs of order six and five, respectively, are given as follows. 

Example 3.8. Consider the cyclic graph on six vertices. Then we have al- 
ready encountered a real regular cyclic sequence in Example l3.3l (0, 4, 2, 3, 1, 2) , 
with corresponding invariant 7 = 10/9. Other complex solutions are indi- 
cated in the figure below. For the left-hand hexagon, we have 7 = 1 which 
corresponds to a 2-colouring of the graph; for the middle one, 7 = 2/3 which 
corresponds to two circuits of a triangle; and for the right-hand one, 7 = —2 
which corresponds to the position function of a regular hexagon. 




1 1 



There is the trivial solution with 7 = 2, when the function (p has value 
on any connected set of either three or two vertices of the hexagon and 1 on 
the connected complementary set of vertices, so for each vertex x, there is at 
most one neighbouring vertex y with (p{y) — (p{x) 7^ (see Section [S]). This 
exhausts all possible values of 7, so the geometric spectrum (the possible val- 
ues of 7) equals the set {—2, 2/3, 1, 10/9, 2} (cf. Section[8|). There is another 
path we have not considered, namely — )• 1 — )• ^ -|- -^i — )• — )• — ^ -|- -^i — )• 
^ + ^1^0, with exterior angles 27r/3, 27r/3, 27r/3, -27r/3, -27r/3, -27r/3, 
respectively. However, this gives the same value 7 = 2/3, corresponding to 
all exterior angles equal to 2tt/3. 

Example 3.9. Now consider a cyclic graph on five vertices. Once more, 
there is the trivial solution with 7 = 2, where we colour the graph with two 
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colours on complementary connected components of three and two vertices. 
There are also two more solutions as indicated in the figure below. 



In these two examples, the angle 9 is given by 27r/5 for the regular pentagon 
on the left, and by 47r/5 for the regular star pentagon on the right. Since 
cos(27r/5) = (^/5-l)/4 and cos(47r/5) = -(\/5 + l)/4, by ([6]), this gives the 
values 7 = — 2/\/5 and 7 = +2/\/5, respectively. These two values together 
with 7 = 2 as above, exhaust the geometric spectrum for the cyclic graph 
on five vertices. 

Consider a planar polygonal chain on A'" vertices {xi, X2, • • • , Xiv} with 
edges all of the same length, such that edges can rotate freely about adjacent 
vertices. In a more general context, where edges have fixed length that are 
not necessarily equal, such an object is sometimes referred to either as a 
planar linkage or as a planar polygonal bar-and-joint framework; they are 
studied notably in robot arm motion planning, see for example [30j . 

A well known problem in the study of planar polygonal chains is to con- 
struct algorithms to either straighten the chain if it is open ended (the 
Carpenter's Rule Problem), or, for a closed chain, to deform it into one that 
is convex. If on allows edge crossings, then this was solved by Sallee in 1973 
|27j . To do this without edge-crossings (for an initial configuration without 
crossings) proved more elusive and was solved by Connelly, Demaine and 
Rote in 2003 0. 

In Section [8] we introduce an energy functional in a more general setting, 
which, for polygonal chains is given by 



When A'^ is even, the absolute minimum of £ is zero which is achieved when 
all exterior angles are zt-Tr, so the polygonal chain is in its most compact form, 
the edges superimposed along an interval. In Appendix [Bl we show how the 
configurations of a polygonal chain on N vertices can be parametrized by 
— 3 parameters and prove that the regular configurations (with exterior 
angle constant) are critical for this functional. The functional £ should 
define a gradient flow on this configuration space whereby a polygonal chain 
evolves into a regular configuration. 




N 



£ = ^(1 + cos6ifc) . 



k=l 
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4. Orthogonal projections of regular polytopes 

In |15) . the authors consider the projection of the vertices of the Pla- 
tonic sohds in onto the complex plane, as well as more general orthog- 
onal projections — >• R^^. As a particular case, they establish that if 
zi, Z2, . . . , zn+1 are the orthogonal projections to C of the vertices of a reg- 
ular simplex in R^, then 

(Zl + • • • + ZN+lf = {N + l)izi^ + ■■■ + . 

In particular, it follows that if we view the 1-skeleton of the simplex as 
a graph, then the function which associates the values zi, . . . , zat+i to the 
corresponding vertices satisfies ([T]) with 7 = N/{N + 1). We now prove a 
corresponding result that shows that, as a graph, the projection of the 1- 
skeleton of a regular polytope satisfies equation ([1]), with 7 constant. We also 
develop some notions that enable us later to define distance and curvature. 

A bipartite graph Ki^n, otherwise known as a star, consists of an internal 
vertex xq connected by edges to n external vertices xi, . . . , Xn] there are no 
other connections. In what follows, we will represent a star embedded in 
M.^ with internal vertex located at the origin, by an (N x n)-matrix 

W = (xil • • • \Xn) , 

whose columns are the components of the external vertices. An invariant of 
the star is a quantity which is invariant under orthogonal transformation of 
the ambient Euclidean space M^. For example the quantity 

-trace (l^Vl^*) = -trace = + • • • ||x„|n 

n n n 

is invariant, where denotes the transpose of W; it corresponds to the 

mean of the squares of the Euclidean lengths of the external vertices. We 

consider a bipartite graph -fCi.n embedded in R^ in the following symmetric 

way. 

Let (2/1, . . . jUn) be standard coordinates for R^ {N > 2); write vectors 
as columns for the purpose of matrix multiplication. Let {ei, . . . , ejv} be the 
canonical basis and write In for the x A^-identity matrix. The internal 
vertex xq is located at the origin, while the external vertices 
situated at distinct points in the hyperplane un = c (constant): 

(7) ^'"(c') (^ = 1'---'")' 

We require further that the (A^ — 1) x n-matrix U = {vi\v2\ - ■ ■ \vn) with 
columns the components v^j of W£ (j = 1, . . . , A^ — 1; £ = 1, . . . , n), satisfies: 

n 

(8) UU' = plN-i, 5^^^£ = 0, 
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for some non-zero constant p (necessarily positive), where denotes the zero 
vector in M''^"^ and C/* denotes the transpose of U . 

Any star, which, up to orthogonal transformation of M^, is embedded in 
this way, we will call a configured star. We shall also say that the vectors 
{iTi, . . . ,Vn} form a configuration in call U the associated configura- 

tion matrix and p the configuration invariant. Provided the star does not 
lie in any proper linear subspace, we say that the star is full. If further, 
llx^ll = r (constant) for £ = 1, . . . , n, we refer to the star as regular of radius 
r. An embedding given by ([7]) and ([8]) is referred to as a standard position 
of the configured star. 

Lemma 4.1. Consider a configured star in (N >2) with internal vertex 
the origin connected to n external vertices {xi, . . . ,Xn} > N). Let W = 
{xi\x2\ ■ ■ ■ \xn) be the N x n-matrix whose columns are the components xgj 
of Xi (j = 1, . . . ,N; i = l,... ,n). Then 

n 

(9) WW* = piN + Tuu\ ^xi = y/n{a+~pju, 

1=1 

where u G is a unit vector, p > and p + a > 0. The quantities n,p,a 
are all invariants of the star; the vector u is normal to the affine plane 
containing xi, . . . , x„. 

Conversely, any matrix W = (xi|x2| ■■■ \xn) satisfying ^ determines a 
configured star with central vertex the origin and external vertices xi, . . . , x„. 

Proof. Consider a configured star in standard position given by ([7]) and dS]). 
Set 



V = 

















Vn \ 


V ^ 


c 







and let A : R — )• R be an orthogonal transformation; set x„ = A 
Then W = (xi|x2| • • • |x„) = AV and 

WW* = AW* A* = piN + a{AeN){AeN)* , 

where 

(10) a = nc^-p. 

Furthermore Yl^=i^t ~ ncAe^, which gives the form Q with u = Ae^- 
The independence of the quantities n, p, a of the orthogonal transformation 
A is clear. 

Conversely, suppose we are given an x n-matrix W = (xi|x2| • • • |x„) 
satisfying ([9]). Let A be an orthogonal transformation such that Au = cm 
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and let V = AW. Write 



V 



Vi 

yiN 



V2 

y2N 



VnN 



Then 



(11) VV' = piN + (TCNeN^ 



and 



E 



ViN 



\/n[a + p) eV , 



so that Ylii = ^-iid y^iv = \/re(cr + p). Furthermore, ([TT]) imphes that 
= p + o". In particular 



n 



But then ^ implies that yiA? = y2N 



VnN = yJ{(T + p)/ 



n. 



□ 



It now follows that the function which assigns the values after projection 
of the vertices of a configured star to the complex plane satisfies ([1]) at the 
internal vertex, independently of the position of the star. 

Corollary 4.2. Let W = {xi\x2\ ■ ■ ■ \xn) define a configured star and let 
P : C be orthgonal projection P(yi, . . . , yjv) = yi + iyAr. Then if 

Z£ = P{xi) = Xii + ixe2, we have 



(12) 



n{a + p) 



where p and a are given by In particular, ^ = a/{a + p) is real and 
depends only on the star invariants. 

Proof. Let u = (ui, . . . ,un) be the unit normal to the plane of the star. 
Then for each j = 1, . . . , A^, we have 



Thus 



^^xtj = \/n{(j + p) ■ 



1=1 



\e=i / 



k,e=i 

n{ap){ui' - + 2\UiU2) 
n{ap){ui + m2f , 
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whereas 

n n 

^ = '^{xei'^ - xn + 2\xaxi2) 
1=1 e=i 

= {WW')ii - {WW%2 + 2\{WW%2 

= (T{ui+\U2f. 

The formula now follows. □ 

In Section [6] we will consider the problem of establishing a converse to 
this corollary. 

Examples of configurations of points vi, . . . ,Vn £ M^^-*^ which satisfy the 
criteria of ^ are as follows. In R, any set of points not all zero distributed 
along the real line with centre of mass the origin form such a configuration. 
In ~ C, we have the n roots of unity: 

(13) = e^"'^/" = 1, . . . n) . 

For n > 3, we have p = ^"^^ cos2(27r£/n) = ^"^^ sin2(27r^/n) = n/2 and 
the coefficient 7 determined by (|12p is given by 

(14) 7 = V- 

This configuration is regular. We note the following elementary fact. 

Lemma 4.3. In M^, every configuration of three points is regular. 

Proof. Consider a configuration of three points in given by Write 
the entries of U as follows: 



U 



Ui U2 U3 
Vi V2 V3 



By a rotation, we can suppose the further normalization: ui = 1, fi = 0. 
This gives the equations; 

U2V2 + nsfs = 
Then this has the solution, unique up to the sign of the square root. 



U 



1 _i _i 

2 2 

n _ 

^2 2 



Since the lengths of the columns of this matrix are all equal to 1, then the 
configuration is regular. □ 
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However, there exist non-regular configurations in M?. Consider the fol- 
lowing set of four points, expressed as the columns of the corresponding 
configuration matrix: 



(15) 



-11-1 1 _ 

y X fi -X -fi J 

where A and /U are real non-zero constants. Then the set forms a configu- 
ration if and only if A^ -|- /x^ = 2. This is regular only when A^ = /i^. A 
non-regular star determined by this configuration will appear later in Ex- 
ample 16.41 In the examples of configurations that follow, we express their 
vertices as the column vectors of the corresponding configuration matrix. 
In M^, the vertices of a tetrahedron placed symmetrically at the points: 

/ -1 -1 1 1 
-1 1-1 1 
\ -1 1 1-1 



(16) 



satisfy conditions ([5]), as do the six vertices of an octahedron: 



and the twelve vertices: 



(17) 



±1 











±1 











±1 





±1 


±A 


±1 


±A 





±A 





±1 



where A is any real constant. In the case when A 



2 ' 



these form the 



vertices of an icosahedron. In M 
vertices satisfying ([8]) given by: 



there is a further configuration of twenty 



(18) 



V 






±A 


±A-i 




±A-i 





±A 




±A 


±A-i 





±1 / 



1+V5 
2 ' 



these form the 



where A is any real constant. In the case when A 
vertices of a dodecahedron (see [13j, §3.8). 

In M^~^, the configuration of 2(A^ — 1) vertices of the cross-polytope: 

/ ±1 ••• \ 
±1 ••• 



(19) 



\ 







±1 



satisfies 
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In order to construct a configuration of vertices of a regular A^-simplex 



in {N > 3) satisfying ([8]), we proceed inductively starting with the 



configuration ()16p in M . 

Suppose we are given vectors vi, . . . ,vn in which label the ver- 

tices of a regular {N — l)-simplex in such a way that, on letting U denote 
the {N — 1) X A'^- matrix whose columns are formed from the components of 
Vi (i = 1, . . . , N), the following conditions are satisfied: 



N 



1=1 



0, UU' = plN-i, 



for constants p and a satisfying p/a'^ = N/{N — 1) (which is the case for 
the 3-simplex of ()16p ). Note that the condition that the vectors do indeed 
correspond to the vertices of a regular (A^ — l)-simplex is that they are of 
equal length and the angle between them has cosine —1/{N — 1). Form the 
following set of iV + 1 vectors in M^: 

/ A^fj \ / cr \ / cr \ 

(20) 



Wo 







Wl 



Vl 



, WN 



VN 



and let W be the matrix whose columns are the components of wi for 
0, 1, . . . , A^. Then it is readily checked that 



N 
1=0 



0, WW^ = piN, 



aN 



VA^2_i 



now with p/ 



aN 



{N + 1)/N, as required. Furthermore, the angle 



between any two Wj and Wk {j ^ k) has cosine —1/N. In particular, the 
vectors correspond to the vertices of a regular A^-simplex and satisfy ([8]). 

To conclude this list of examples of configurations, we note that if (i7i|'V2| ■ ■ ■ 
is a configuration in M^~^ with invariant p, then both 



/ Vl 


V2 




Vn 


' ) 


\ ^ 


c 




C 


—nc J 



where c = p/ \J n(l + n), and 



/ ^1 


V2 




Vn 







V 










V 2 


4) 



are configurations in with the same invariant p. The first of these gen- 
eralizes the inductive construction of the vertex configuration of the regular 
simplex given by ([20]) . 

Before proceeding, we assemble the information we require concerning 
regular polytopes, which can be found in the classical text of Coxeter jl3| . 
The regular polytopes have a symbolic representation in terms of the Schlafli 
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symbol. This can be defined inductively as follows. For a regular polygon 
with p edges, one assigns the symbol {p}. A regular star polygon which 
winds m times around its center is denoted by the fractional value {p/m}, 
so for example, the second graph of Example 13.91 has Schlafli symbol {5/2}. 
A regular polyhedron whose faces are polygons of type {p} which has q such 
faces joining around a vertex has symbol {p,q}- A regular 4-polytope (or 
polychoron) with highest dimensional cells polyhedra of type {p, g} having 
r such cells joining around an edge has symbol {p, <?, r}, and so on. 

A regular polytope in MJ^ with Schlafli symbol {p, g, . . . , s, t} is charac- 
terized by its highest dimensional cells which are regular polytopes with 
symbol {p, ...,s}, and its vertex figure^ which is a regular polytope of 
type {q,r ■ ■ ■ ,i} in M^^-*^ obtained by fixing a vertex xq and constructing 
a polytope in an affine {N — l)-plane, whose vertices are points half way 
along each edge emanating from xq. For example, the dodecahedron has 
Schlafli symbol {5,3}; it is made up of three pentagonal faces around each 
vertex; its vertex figure is the triangle obtained by moving half way along 
each edge emanating from a particular vertex and joining these points by 
edges which traverse each of the three corresponding faces. The 600 cell is a 
4-dimensional regular polytope with 120 vertices, 720 edges, 1200 faces and 
600 tetrahedral cells; it has 5 tetrahedra joining around each of its edges 
and so has Schlafli symbol {3, 3, 5}. Its vertex figure is an icosahedron with 
symbol {3, 5}. 

In two dimensions, the regular polytopes are just the regular polygons and 
star polygons, as discussed in Section [3] (see Example [32] for an illustration of 
a star polygon on five vertices). A rich variety of regular polytopes exists in 
three and four dimensions. In three dimensions, there are five convex regular 
polytopes (polyhedra): the tetrahedron, octahedron, cube, icosahedron and 
dodecahedron, often known as the Platonic solids; in addition there are four 
non-convex regular polytopes called star polyhedra, with Schlafli symbols 
{|,5}, {5, |}, {|,3}, {3, |}. 

The star polyhedra can be constructed from the convex polyhedra by a 
process known as stellating and faceting. For a polygon, stellating consists 
of maintaining the edges and extending them until they connect in new 
vertices; faceting on the other hands consists of maintaining the vertices and 
inserting new edges in an appropriate way. The procedure for polyhedra is 
similar: either the faces are extended to create new vertices or the vertices 
are maintained and new faces are constructed (see ^13j, Chapter 6). In 
particular, the ensemble of positions of the vertices of a star polyhedron is 
always congruent to those of a convex polyhedron. 
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In dimension four, there are sixteen regular polytopes, six of them convex. 
An important property to note is that it is precisely the three dimensional 
regular polyhedra that arise as vertex figures of the four dimensional regular 
polytopes. 

In dimensions five and above, there are just three kinds of regular poly- 
tope, the regular simplex with Schlafli symbol {3,3, . . . ,3}, the cross poly- 
tope with Schlafli symbol {3, 3, ... , 3, 4} and the measure polytope, or hy- 
percube with Schlafli symbol {4, 3, . . . , 3}. 

The above discussion leads to the following consequence. 

Theorem 4.4. Let V = iV^E) he the graph given by the 1- skeleton of a 
regular polytope and let P : C be any orthogonal projection of the 

ambient Euclidean space. Let ip = P\v : ^ — )■ C. Then ip satisfies ^ with 
7 constant. 

Proof. It suffices to note, that in each case the vertex figures with respect 
to a particular vertex xq can be positioned in such a way that the vectors 
of the corresponding star satisfy the criteria of Corollary 14.21 In the case of 
a polyhedron, the vertex figure is a polygon whose vertices can be placed as 
in (|13|) . For a 4-dimensional regular polytope, the vertex figure is a regular 
polyhedron and, as the case by case list above shows, their vertices can also 
be placed in the required way. Similarly, for the cross-polytope, whose vertex 
figure is another cross-polytope of one dimension lower; this is dealt with 
by the vertex placement (jl9|) . The hypercube and the regular A^-simplex 
both have vertex figure a regular — 1-simplex, whose required placement 
is given by the inductive construction of the vectors given in ()20p (with N 
replaced by — 1). □ 

5. Invariant frameworks 

A framework T in is a finite collection of points {xi, . . . , Xn\ which 
correspond to the vertices of a graph whose edges are straight line segments 
joining the vertices. We shall view a framework as a graph and continue to 
use the notation T = {V, E) to distinguish the vertices and edges. The edges 
are often called bars and such bar frameworks have been much studied in 
respect of questions about rigidity |10j . Thus one is interested in whether it 
is possible to deform the bar framework whilst preserving the length of the 
edges. If the only such deformations correspond to rigid motions, then the 
bar framework is called rigid. The planar linkages of Section [3] (other than 
the triangle) provide examples of non-rigid bar frameworks. 

We have seen in the last section that the framework corresponding to the 
1-skeleton of a regular polytope has the property that the function which, af- 
ter an orthogonal projection — )• C associates the corresponding complex 



A CLASS OF QUADRATIC DIFFERENCE EQUATIONS ON A FINITE GRAPH 23 



values to each vertex, satisfies ([T]). Furthermore, it does so in an invariant 
fashion with respect to orthogonal transformation of the ambient Euclidean 
space. In this section, we wish to determine more general frameworks for 
which these properties hold. Amongst the examples are ones that are "flex- 
ible" in the sense that they can be deformed continuously through solutions 
to ([1]) without disturbing the value of 7. 

Let J- = (V, E) be a framework in and let P : — )■ C be some 
orthogonal projection of the ambient Euclidean space. Let ip = P\y : V ^ 
C. We say that the framework is invariant if (i) 93 satisfies ([1]) for some 
7 : y — )• M; (ii) if A : — )• is any orthogonal transformation, then 
if A = P o A\y : V ^ C also satisfies ([TJ with corresponding 7^ satisfying 
'jAiAx) = 7(2;) for all x €z V. Note that if (i) is satisfied, this always remains 
the case if the framework is translated or dilated, since this just corresponds 
to a normalization ip 1— )• Xip + (now with A real). Then condition (ii) 
implies invariance by a transformation x 1— t- XAx + b, where A > is real, A 
is any orthogonal transformation and b £ is any vector. 

By Theorem 14.41 any framework given by the 1-skeleton of a regular 
polytope is invariant. But do there exist other invariant frameworks that 
arise as the 1-skeleton of non-regular polytopes? It turns out that the answer 
to this question is yes. In Section [UJ when we consider connected graphs 
which satisfy ([1]) as components of a more complex system, we will view such 
objects as "physically" significant. We now construct families of examples. 
In order to do this, we generalize the configured stars of the last section. 

Consider a star in M.^ with internal vertex the origin and external ver- 
tices located at the points xi,X2, ■ ■ ■ ,Xn- As in the last section, it is often 
convenient to represent the star by an x n -matrix whose columns are 
the components of the vectors Xi (i = 1, . . . ,n): 

W = {Xi\x2\ ■ ■ ■ \Xn) ■ 

Consider the projection P : — )• C given by P{yi, ■ ■ ■ , un) = yi + 'W2- Let 
A : —7- be an orthogonal transformation and set Z£ = P o A{x£). We 
say that the star is invariant if it satisfies the equation: 

(n \ 2 n 
£=1 / e=i 

with 7 real and independent of the transformation A. By Corollary 14.21 ^ 
configured star is invariant. However, there are more general invariant stars. 
One family is given by the following lemma. 



24 PAUL BAIRD 

Lemma 5.1. The star in "M? with 2r external vertices represented by the 
columns of the 3 x (2r) -matrix 



Xl 


X2 ■ ■ 




Xl 


X2 ■ ■ 


Xf \ 


Si 


S2 ■ ■ 




-Si 


-S2 ■ ■ 


• -Sr 


h 


t2 ■■ 




-h 


-t2 ■■ 


■ -tr / 



where the vectors s = (si, . . . , Sr) and t = (ti, . . . ,tr) are orthogonal and of 
the same length, is invariant. 

We omit the proof, which runs along similar lines to that of the theorem 
which follows. 

The lemma enables us to construct a family of invariant polytopes in M^, 
as illustrated in the following figure. 



We take a regular polygon on n vertices in the yiy2-piaiie with centre 
the origin and join each of these vertices to two symmetrically placed points 
(0, 0, ±6) along the ys axis, where 5 is a constant to be determined. In the 
first instance, suppose the vertices of the polygon are located at the points 
^2km/n ^ _ Q^]^^__ ,n— 1. Fix attention on the star centred on the vertex 
at (1,0,0) G M^. After translating to the origin, the matrix of the star is 
given by 

/-I -1 + cos^ -1 -1 + cos^ 
sin^ -sin^ 

n n 

\ h -h 

We then see that the uniquely determined choice b = sin ^ (up to sign) 
guarantees that this matrix has the form of Lemma 15.11 We may calculate 
7iat at this lateral vertex, to give: 

2(l-2cosf + 2cos^^) 
^^^^ (2-cosf)2 • 

On the other hand, the stars at the apexes, are regular configured stars 
which therefore satisfy (j2T|) . with 7apcx given by (dH): 

2 sin^ ^ - 1 

(23) 7apex = 2sin2 2^ ■ 
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The two values of 7 coincide precisely when n = 4, in which case 7 = 1/2; 
then the figure is regular and corresponds to the cross-polytope. It is readily 
checked that any other choice of b destroys invariance. We now generalize 
this construction. 

Consider a framework F in given by the 1-skeleton of a regular poly- 
tope. We suppose the polytope is centred on the origin in R^. Define the 
double cone on F of height b to be the new framework in M^"^^ obtained 
by adding two vertices symmetrically placed at the points (0, 0, ... , 0, ±6) S 
M^"^^ and adding edges joining each of these points to each vertex of the 
polytope. The example discussed above represents a double cone on a regu- 
lar polygon. (Single) cones on frameworks have been considered in respect 
of rigidity problems by Connelly and Whitely [11] . 

Theorem 5.2. Given the framework F corresponding to the 1-skeleton of a 
regular polytope, then there is a unique double cone on F which is invariant. 

In particular, if the regular polytope is convex, then the double cone 
corresponds to the vertex and edge set of a new convex polytope of one 
dimension greater. In general this is only regular if the original polytope is 
the cross-polytope. Indeed, the degree is no longer constant in general. For 
example, the double cone on the dodecahedron has two vertices of degree 
20 and twenty vertices of degree 5. These polytopes generalize the regular 
polytopes, as being invariant by rigid transformation in M^+i in respect of 
equation ([T|). 

Proof of theorem : Consider a regular polytope in M^. Fix attention on one 
particular vertex x which we suppose located at the origin in R^. We may 
suppose the coordinates chosen so that the vertex figure at x is given by the 
configuration U = {vi\v2\ • ■ ■ \ vn) M^~^, where M^"^ is embedded in 
as the first — 1 coordinates. In particular, the star on the vertex x is given 
by the matrix: 



for some non-zero constant c. 

Now the line from the origin to the centre of the configuration passes 




through the centre of the polytope, which is at some point 



a 







) 



where 



G M^-i 



is the centre of the configuration given by U. 
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In order to construct the double cone, we situate the polytope in the plane 

/ 

Un+i = in and add two more vertices at the points 



. These 



\±b 



two vertices are to be added as new external vertices to the star centred on 
G M^+\ whereby the star matrix now becomes: 



(24) 



S 





V2 




Vn 







c 


C 




C 


a 


i) 


\ 










b 





This is not in general a configured star. We wish to confirm its invari- 
ance under orthogonal transformation with respect to ()2ip (with zi, . . . ,Zn 
replaced hy zi, . . . , Zn, Zn+i, Zn+2 as defined below). 

Let A : M^"*"^ — )• M^"'"-'^ be an orthogonal transformation and as usual 
set A = {ars)i where now r, s G {1, 2, . . . , + 1}. Let j^k range over 
the indices {!,... ,A^ — 1}, so that the components of vi are given by vij. 
We compute the first two rows of AS to determined the complex numbers 
Z£, Zn+i, ZnJ^2) where 1 < ^ < n, after projection P : — t- C given by 
P{yi,. . . ,yN+i) =yi+W2- 

N-l 

c{aiN + ia2Ar) + ^ {aij + 'm2j)vij 
i=i 

a{aiN + ia2N) + b{ai^N+i + ia2,Ar+i) 
a{aiN + ia2N) - b{ai^N+i + ia2,Ar+i) 



Zn+2 



On recalling that for each j = 1, . . . , N — 1, the sum V£j = 0, we have 

n 

(25) Zn+i + Zn+2 + '^ze = {nc + 2a){aiN + ia2N) ■ 
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Now 

N-1 

zi' = c^iaiN + ia2Nf + 2c{aiN + ia2N) ^ (aii + '^a2j)vij 

j=i 

N-l 

+ + ^o-2jfvij'^ + 2^(aij + ia2j)(aife + ia2k)vtjVik ■ 

j=l j<k 

But since U = {vi\v2\ • • • \vn) is a configuration, from ([8]), we have 

n 
£=1 

for all j, /c = 1, . . . , — 1. It now follows that 

n 

On comparing with (j'25|) . we see that (j21|) is satisfied with no dependence 
on ^ if and only if 6 = p/2. In this case 

_ (n + 2)(2a2 + nc2 -/)) 
~ (nc + 2a)2 ' 

where 7iat refers to the lateral value of 7, namely the value at one of the 
vertices of the regular polytope. 

For the two vertices corresponding to the apexes of the double cone, the 
invariance if given by Corollary 14.21 Indeed, as observed in the last section, 
the vertices of any regular polytope V, convex or not, form a configuration. 
If we let pp denote the invariant of this configuration and let m denote the 
its cardinality, then from (|10p . we have 



a = mh - p-p = — - p-p . 

It follows that 

/ORN mp - 2pp 

(26) 7apcx = — : = ■ 

a + pp mp 

This completes the proof of the theorem. □ 

We have computed the two values of 7 in the above proof so as to be 
able to apply them to examples. Note that we are at liberty to normalize 
the parameters o, c, p as we wish. For example, for the double cone on a 
regular polygon of n sides, c = 2asin2 ^ and p = 2a} sin^ This confirms 
the values of 7 computed in the example preceeding the theorem. 

If our initial polytope is the cross-polytope in R-'^, then we have m = 2N 
and n = 2{N — 1). If we normalize so that a = 1, say, then c = 1 and 
p = Pp = 2. It follows that the unique value of b which gives an invariant 
double cone is 6 = a/ p/2 = 1, as required, since the double cone on the 
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cross-polytope is another cross-polytope of one dimension higher. One easily 
checks that 7apex = 7iat in this case. 

Before moving on, we estabhsh one more instance of invariant frameworks 
which has imphcations for questions of "rigidity" . 

Proposition 5.3. Let {vi\v2\ ■ ■ ■ \vn) be a configuration in M^. Consider 
the framework T given by the complete graph on this configuration. Then 
there is a unique double cone on T which is invariant. 

Proof. We use the same notation as in the proof of the above theorem, except 
that now each vector V£ lies in rather than in M^-\ so the indices j, k 
range over 1, . . . , A^. The invariance at the apexes is guaranteed by Corollary 
14.21 We therefore fix attention on one of the lateral vertices, say vi. Perform 
a translation so that vi is located at the origin. If we suppose the height of 
the double cone is given by 6, then the star matrix at vi is given by 



/ V21 

V22 



vn 

Vl2 



V31 
V32 



Vn 

Vl2 



V2N - VlN V3N - ViN 

V 



Vnl 
Vn2 



Vn 
Vl2 



VnN - VlN 





-Vn 

-Vl2 



-VlN 

b 



-Vn \ 

-Vl2 
-VlN 

-bj 



Now consider the effect of an arbitrary orthogonal transformation A = (oij) : 
R-'^"'"^ — )• M^"*"^ on the star followed by projection to the first two coordinates. 
This gives the corresponding complex numbers: 

TV 

Z£ = ^(aij + ia2j){vij - vy) {£ = 2,... ,n) 



N 



b{"-i,N+i + ia2,Af+i) - + ^a2j)vij 



Zn+2 



-bi<^i,N+i + ia2,Ar+i) 



TV 

^{aij + ia2j)vij 
i=i 

But since the vectors vi form a configuration, we have Yll=2 ^^i 
allj = 1,...,A^. Thus 

n TV 
Zn+1 + Zn+2 + '^^1 = -{u + 2) ^(a^ + ia2j)vij . 



-vij for 



For the sum of squares, we obtain 

2 I ^ 2 , Y^n 2 
Zn+l + Zn+2 + 2-^£=2 



+ '^a2jf{vij^ - b'^) + 4V<fc(«ii + i«2i)(aifc ia2kVij)vik 



+ E"=2 Ej=i(«ii + i«2i) {vej - vij) 

+2 I]"=2 Yyj<ki"-ij + ia2j)(aiA: + ia2k){vej - vij){vek - vik) 
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If we let p denote the configuration invariant, so that ^^^i vijvik = P^jk^ 
then this gives 

n N 

Zn+i' + + y^^zi^ = ^(oij + ia2j)^ [(n + 2)ui/ + - 26^] 

1=2 j=l 

+2(n + 2) ^(aij + ia2j){aik + ^a2k)vijVik ■ 
j<k 

On the other hand 

^ Zn+l + Zn+2 + X]"=2 = 

+ 2)2 ^Y.f^^iaij + ia2j)2wi/ + 2X;j.^fc(aij + ia2j)(aifc + ia2fc)'t^ii^'ifc} • 

It now follows that equation (|2ip is satisfied (with zi , . . . , 2;„ replaced by 
Z2, • • • , 2;„+2) with 7 independent of j4 if and only if 26^ = p, in 

which case: 



(^Zn+l + Zn+2 + = ^n+l^ + Zn+2^ + 

This completes the proof. 



n + _ 

1=2 1=2 



n 



The significance of this proposition is that it provides examples of frame- 
works that are "flexible" in respect of equation ([T|), by which we mean the 
edge lengths may change, but the framework can be deformed continuously 
through solutions to ([1]) while maintaining invariance. Configurations for 
which this is the case are given by the examples (fT5]l , (fT7|) and (fT8]l of Sec- 
tion m In each case, there is a 1-parameter family of configurations which 
yield invariant double cones. 

The last formula of the proof shows that the lateral value of 7 is given by 

n-M 

7lat = -T, ■ 

n + 2 

At an apex, the value of 7 is obtained from ()10p . In fact 

a = nh'^ - p = {n- 2)b^ , 

so that 

a n — 2 

Tapex — ; — • 

a + p n 

We note also that 7iat and 7apex depend only on the cardinality n of the 
configuration. Although these two values are never equal, remarkably, when 
vi, . . . ,Vn are the vertices of a regular polygon in the plane, the addition of 
one more edge joining the two apexes yields an invariant framework in M'^ 
with 7 = (n -|- l)/(n -|- 2) constant. The addition of this extra edge makes 
the framework into a complete graph on n -|- 2 vertices. 
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Say that a framework is immersed if all vertices are distinct; say that it 
is embedded if it is immersed and no two edges intersect except at their end- 
points. Then the complete graph on n + 2 vertices can always be embedded 
as an invariant framework in M"+^, its vertices and edges corresponding to 
the 1-skeleton of a regular n + 1-simplex. 



Corollary 5.4. There exists an invariant immersed framework in with 7 
constant corresponding to the complete graph onn+2 (n>2) vertices. There 
exists an invariant embedded framework in M" with 7 constant corresponding 
to the complete graph on n + 2 (n >2) vertices. 

Proof. Let vi, . . . ,Vnhe unit length vertices of a regular polygon in the plane 
and construct the unique invariant double cone of the above proposition. 
Since p = n/2 (Section [3]) , its height is given by 6 = -y/n/^- Lemma [A. II 
of Appendix |A] shows that the addition of another vertex at a distance x 
(7^ —nb) along the axis of symmetry of the cone connected to the apex, 
yields another invariant framework, with new 7 at the apex given by (|47|) . 
that is by 

_ (n+ +n5^7apex) 

Tapex / , xn ) 

[x + nby 

where 7apex = ~ 2) /n is the value of 7 prior to the addition of the new 
vertex. But now substitution of the values b = -y/n/2 and x = \/n yields 
precisely the value 7apex = (f^ + + 2), which is the value at a lateral 

vertex. Furthermore, the new edges connects the two apexes. 

For the second part of the corollary, when n > 3 one proceeds as above 
but now with v\^. . . ^Vn the vertices of a regular (n — l)~simplex in M*^"^. 
Indeed, as noted after the above proposition, the value of 7iat depends only 
on the cardinality of the configuration. The resulting double cone is clearly 
embedded, since when n > 3, the centre of mass of the regular simplex is not 
contained in any of its edges. But the edge joining the two apexes intersects 
the subspace M"~^ precisely at this point. When n = 2, the framework can 
be taken to be the projection of the 1-skeleton of a tetrahedron as indicated 
in the right-hand figure below. □ 

In Appendix |X] we show how to build more intricate frameworks, called 
invariant structures, from the elementary ones constructed above. Examples 
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with 7 constant such as those given by the above corollary, are particularly 
significant in respect of our model of an elementary universe developed in 
Section [9l where we view them as stable particles. 

Given a solution to ([T]) on a connected graph T, the global embedding 
problem is to realise the graph as an invariant embedded framework J-" = 
{V, E) in a Euclidean space in such a way that '^{x) = P{x) for all 
X £ V. The above corollary shows that the solution corresponding to the 
projection of the framework of a regular n-simplex (in R") can be embedded 
in M"~^ for n > 3. Clearly, if T can be realised as a planar graph with non- 
intersecting edges in such a way that ip is the position function associated 
to the vertices, then the solution {T,ip) is globally embedded in M^. This is 
the case for the left-hand solution below, for which 7 = 1/3 (see Example 



In fact, by direct calculation, one can show that this solution admits no 
global embedding in and by the lemma below, it cannot admit any global 
embedding in M.^ for > 3, hence M? is the unique Euclidean space into 
which it embeds. 

Lemma 5.5. Let ip be a solution to ^ on a connected graph T = {V, E) 
such that 7(x) < 1 for all x £ V. Suppose that uq is the smallest degree 
of the vertices of T . Then the solution cannot be realised as an invariant 
embedded framework in any with N > uq. 

Proof. Suppose on the contrary that there exists an invariant embedding in 
with N > Uq. Let x be a vertex of degree no and consider the affine 
subspace S spanned by the edges emanating from x. Let P : — C be 
an orthogonal projection. Then since dim 5 < N, there exists a Euclidean 
motion A : — )• such that P{A{S)) is entirely real. In particular each 
neighbour y oi x has P{A{y)) real. But by Lemma |3.6^ this contradicts the 
hypothesis that 7 < 1. □ 

The right-hand figure above shows a projection of the tetrahedron. This 
can be realised either as an embedded framework in R^, or in M^. On the 
other hand, a solution corresponding to the projection of a cube cannot be 
embedded in the plane, since edges must cross; is the unique space into 



[631): 




—I 
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which it embeds. Corollary 15.41 shows that the complete graph on five ver- 
tices admits solutions to ([1]) with 7 = 4/5 which correspond to embeddings 
in and M^, and only into these spaces. The star graph of Example 13.91 
shows that there exist examples with 7 < 1 which cannot be embedded 
into any Euclidean space: since the degree of each vertex is 2, by the above 
lemma it can only embed into M?, but then edges must cross. In the next 
section, we address the local embedding problem, or lifting problem and 
show that, for a given solution to ([T]), we can always embed a vertex and its 
neighbours in an invariant way. 

6. The lifting problem 

Given a solution to ([1]), at each vertex x, our aim is to construct a 
configured star in some Euclidean space whose external vertices project 
to the points ip{y) — ip{x) {y ~ x) of the complex plane. To do this, we 
establish a converse to Corollary 14.21 We shall refer to the problem of 
constructing such a star as the lifting problem. At a vertex of degree three 
with if holomorphic, this is the Theorem of Axonometry of Gauss |17j . It 
turns out that provided 7 < 1 (n = vertex degree), the lifting problem can 
always be solved. Remarkably, when = 3, the configured star is unique 
up to a sign ambiguity and arises as the minimizer of a natural functional 
determined by ip. Relative edge-length and so relative distance on the graph, 
can now be defined in terms of virtual configured stars. 

Suppose we are given zi, . . . , z„ G C (n > 2) not all zero satisfying 

(n \ 2 n 
^zA =^z/ (7GM). 
e=i / 1=1 

For a given N with 2 < A < n, we wish to construct a configured star 
W = (xi|^2| ■■■ \xn) in with Z£ the orthogonal projection of xi. For 
convenience, write zi = xn + 1x^2 = o-i + i/3^, so that 





1 OLl 


a2 




Oln \ 










Pn 


w = 




X22, 




Xn3 




\ XiN 


X2N 




XnN j 



For A > 3, we are required to solve the system: 

n 

(28) WW* = piN + (Tuu\ '^xe = y^n{a + p) u , 

1=1 

for X£j {I = 1, . . . , n; j = 3, . . . , A), p > 0, cr such that p -1- cr > and 
u E unit, with ^ = a/{a + p). If we require the star to be regular, then 
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we have the further condition: 



Xi\ \ =r 



\/i = 1, . . . , n, for some r > . 



If = 2, then W is determined and for n > 3, the system ()28p is not in 
general satisfied; if = n = 2 it is always satisfied. 

Note also the normalizing freedom; namely, (|27|) is invariant by the re- 
placement Z£ I— Xz£ (A G C). We will discuss the effect of this on the 
parameters of the problem below. 

Theorem 6.1. Given a non-zero solution {zi, . . . , Zn',j} to ()27p with n > 3 
satisfying 7 < 1; then for each 3 < N < n, there is a full configured star 
W = (xi|x2| • • • \xn) in such that zi is the orthogonal projection of xi . 
There are two cases: 

(i) If N = 3 there is a 3 x n-real matrix A and a vector 6 G M^, both 
of which depend only on zi, . . . , Zn, such that the vertices of the star are 
determined by the solutions t £ to the equation At = b; provided the zg 
do not satisfy the identity 



then the star corresponds to the minimal Euclidean norm solution to this 
equation and is unique up to a sign ambiguity. 

(ii) If N > 3, the star is no longer unique, and although it corresponds 
to a solution X to a more general system AX = B, where X is a variable 
n X (N — 2) -matrix and B is an 3 x (N — 2) -matrix which depends only 
on zi, . . . , Zn, it is not in general minimizing amongst solutions. 

In the case when = 3, except for special cases, the theorem says that the 
configured star arises as the absolute minimum of a functional determined 
by the solution to ([27]) . The special solutions which satisfy (f2U|) are char- 
acterized in Lemma 16.31 below. Note that we do not insist that Z£ 7^ 0, 
which is necessary to determine 7; in fact in this case we have a family 
of stars for different choices of 7 subject to 7 < 1. In order to prove the 
theorem, we first of all establish some preliminary identities. 

Lemma 6.2. // the system ([28]) is satisfied, then we have the following 
equalities: 



(29) 




e 



e 



(30) 




(31) 
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In particular, we have the identity: 

e e e 

Furthermore, provided 7 < 1, then p > 0. 

Proof. Identity (j30p follows from the second equation of (j28p . Then the first 
equation of ([28]) implies that 

which gives the first identity of (j31|) . The second identity follows similarly. 
The inequality p > is an immediate consequence of (l5|). □ 

These are necessary conditions that must be satisfied by any configured 
star that projects to a solution of ()27p . 

/A \ 



Write vectors in column form and set a = 



and 



\ On J 

1 E M"' to be the vector having 1 in each of its entries. 

In the context of equation ()27p . by normalization we mean the freedom 
to multiply each zi by a non-zero constant A G C; this leaves the equation 
invariant. So we no longer allow the more general freedom ip 1— )• Xip + fj, asso- 
ciated to equation ([1]) since that disturbs the requirement that (^z-^, . . . , -^n) 
represents a solution which is zero at the internal vertex of a star. The 
following lemma concerns normalizations that are useful in the proof of the 
theorem. 

Lemma 6.3. (i) If Z£ ^ 0, then there is a normalization with respect to 
which Zi = 1, in which case: 

(32) J^a£ = l, X1'^^ = 0' Z]«^'^^ = °- 

e e e 

(ii) For a non-trivial solution to ()27p with 7 < 1, suppose there exists a 
linear combination of the form aa + bf3 + cl = for a, 6, c G M not all zero. 
Then there exists a normalization factor A such that Xzi = ai + if3, with 
a£ gM (i = 1, . . . ,n) and with /? G M constant. Furthermore /3 7^ and the 
identities: 

(33) ^a£ = and 7 = 1 " ^ X] ' 



are satisfied. In particular (j29p holds. 
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Conversely, if (j29p holds, then there is a linear combination of the form 
aa + 6/3 + cl = for a, 6, c G M not all zero. 

Proof, (i) This is trivial; if Yle Z£ = fi ^ 0, then we multiply each Z£ by 
A = fJ-^^. Then after normalization 

n n 

so that J2e <^il^i = 0- 

(ii) Suppose there exists a linear combination of the form aa + 6/3 + cl = 
for a, 6, c G M not all zero. If 6 7^ 0, then 

z,= (l-i^) ae-il 

so that 

b{b + ia)zi = {(a^ + 6^)a£ + ac} — icb , 

which is of the required form. If on the other hand 6 = 0, then since 7 < 1, 
by Lemma 13.61 a 7^ and —iz£ = f3i + i(c/a), which is once more of the 
required form. 

Suppose now that Z£ has the form Z£ = a^ + i/3 with /3 G M constant. First 
note that /3 7^ 0, otherwise we would contradict Lemma 13.61 once more. The 
identities ()33p now follow by taking the real and imaginary part of ()27p . It 
is straightforward to check that these imply (|29p . 

Conversely, suppose that ([29]) holds. If Yle-^^ ~ 0' then from ([29]) . 
^^jz^P = which implies each Z£ vanishes, contradicting our hypothesis. 
Hence zg / 0. Now normalize so that zi = i. Then 

n n 

and 

n 

7 = -nJ^(a,2_^/). 

But now (j29l) implies that 

n n ( ^ \^ 

nJ^/3,2-l = o ^ n^/3/-K]/3, =0, 
i=\ £=1 \e=i J 



which, by the Cauchy-Schwarz inequality implies that I3i = (3 constant. 
Since 7 < 1, we must have /3 7^ and there exists a non-trivial linear 
combination: /? — /31 = 0. □ 
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Proof of Theorem 16.11 Let {zi, . . . ,z„;7} be a non-trivial solution to ([2 
satisfying 7 < 1. Set 



(34) 
and 

(35) a = 

Then from (1271) we have 



^7 + p = p/(l - 7) > 0) 



2 7 



and 



E«^)(E/30 = E"^/3^' 



so that the identities (ISTTl hold: 



i e I e 

Define ui and U2 by ([30]) . We are required to solve ([28]) for xij {I = 
1, . . . , n; j = 3, . . . , A^) and for U3, . . . , un satisfying ui^ + + + • • • + 
tiTV^ = 1- It is convenient to rewrite (|28|) as follows. 
Set 



A 



Then 



V 1 1 



Or 



/3n , ^ 



/ 2;i3 Xi4 
2^23 X24: 



XlN \ 

X2N 



AX 



/ (a,X3) (a,X4) 



{a,XN) 



where {d,Xj) denotes the Euclidean inner product of the vectors a = 

, for j = 3, . . . , N (similarly for (/3, Xj)). Then 









( ^ij 


\ 






and Xj = 






\ On 


) 




\ Xnj 


J 
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solving (j28p is equivalent to solving 



AX = B : 



(36) 



auius 



auiUN 

CrU2UN 



\ ^/n{a + p) Us ^Jn[a + p) U4 ■■■ ^Jnia + p) un 



(JU2 



Us U4 



subject to the constraint: 

/ {Xs,Xs) {Xs,X4) 



(37) X^X 



{X^jXs) (X^jXi) 



\ (XnjXs) {Xn,X4) 
p + (JUs^ ausUi 



(JU4US p + cn4 

y auNUs cruNU4^ 
( us \ 



pi IS! -2 + (y 



U4 



{us U4 



UN 



{Xs,Xn) \ 

{X4, X]\f) 

{Xn,Xn) J 

ausUN \ 
au^UN 

P + gun'^ j 
■ un) 



\ UN J 

It is important to note the sign ambiguity: the equations are invariant under 
the simultaneous replacement of Uj by —Uj for j = 3, . . . , A'" and of X by 
—X. This ambiguity represents two choices for the configured star. In order 
to study solutions to p6]l , we will employ the Moore-Penrose pseudo-inverse 
of the matrix A [23j. 

In general, if ^ is an m x n -matrix with m < n such that the rows of 
A are linearly independent, then AA^ is invertible and the Moore-Penrose 
pseudo-inverse is the n x m -matrix 

A+ := A\AAY^ . 

If is an m X p -matrix and X a variable n x p -matrix, then Z = A^B is 
a solution to the equation 

AX = B; 

furthermore, it satisfies < ll^ll-F fo^^ solutions X, where 



\X\ 



i=i 3=1 



trace (X*X) , 
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is the square of the Frobenius norm of X. The general solution is given by 

X = A^B + {In-A+A)V, 

where V is an n x p -matrix which can take on arbitrary values. We now 
return to the system (j36p . 

Suppose first that the rows of A are linearly independent. Set Z = A^B. 
Then AZ = B. In order to study the constraint (|37p . we calculate Z^Z: 

Z^Z = {A+BfA+B = B\A+YA+B 
= B\A\AAYyA\AA^)-'^B 
= B\{AAYyAA\AA^)-^B 



B\AA 



B . 



Then from (1301) and (I3TI1. 



AA' 



auiU2 



auiU2 
p + (Tn2^ 



\ y^n{a + p) ui ^Jnia + p) U2 
with determinant: 

(38) det(^^*) = np^(y\ - - 

The inverse is given by: 



\Jn[o + p) ui 
\Jn{a + p) U2 
n 



U2 



{AA 



1 



np{l — ui' 



U2^ 



n(l - n2^) 
nuiU2 



\ -y/n(a + p)ui 
We then compute to find 



-y^n(a + p)ui 
-y/n{a + p) U2 
■\/n{a + p) U2 /3 + fT(ni^ + ti2^) 



nuiU2 
n(l — ui^) 



(39) 



P 



/ us \ 

U4 



{U3 Ui 



UN 



1 — Ui^ — U2 

\ UN J 

On comparing (j39p with ()37p . we see that we always have equality of the 
right-hand sides when = 3. In particular, the solution X = Z which 
satisfies the constraint minimizes ||^||f among solutions X to AX = B and 
so is unique up to the sign ambiguity referred to above. (If also n = 3, then 
there is one and only one solution to AX = B; this solution automatically 
satisfies the constraint). On the other hand, for > 3, the solution X = Z 
never satisfies the constraint and we are required to look at a more general 
set of solutions. 

If an n X — 2) -matrix Y satisfies the homogeneous equation AY = 0, 
where is the 3 x (N — 2) zero matrix, then X = Z + Y solves the equation 
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(j36p . We now look at the homogeneous equation. Write the column vectors 
of y as Yj (j = 3,...,iV). 

Suppose first that Yli 0- Then by Lemma 16.31 there is a normal- 
ization with respect to which the conditions (j32|) are satisfied. Note that 
in this normalization U2 = and ui = 1/ \ln(a + p) = ^J^jna . Consider 
R"' with coordinates (ti, . . . ,tn)- Let 11 be the (n — l)-dimensional hyper- 
plane + • • • -|- tn = with inner product induced from the standard one 
on M"'. Then /3 G 11. Let vr : M" — )• LI be orthogonal projection and consider 
7r(a); this vector (which could vanish) is orthogonal to /?. Then a necessary 
and sufficient condition that Y satisfy AY = 0, is that the column vectors 
Yj (j = 3, . . . , A^) lie in LI and are orthogonal to both /? and 7r(a). Let us 
now consider the constraint (j37|) for X = Z + Y. 

Then 

x^x = z^z + y*z + Z^Y + Y^Y . 



Given (|39p . we see that X satisfies (j37p if and only if 



(40) plN-2 = 2^ \ 2 

W3 -I 1- UN 



( ns \ 

\ UN J 



{u3Ui ■■■ un)+Y'Z+ZY+YY 



Now a calculation shows that 



(41) Z = A+B 



n(l — ui' 



W2 



02 (^2 1 



-nui 
-nu2 



U3 



UN 



Recall that in we use coordinates (yi, . . . ,yN)- By a rotation of the 
(A^— 2)-dimensional subspace spanned by (ys, . . . , t/at), we may suppose that 
the vector {u^, . . . , un) is directed along d/dy^; thus = = ■ ■ ■ = un = 
and u^^ = 1 — ui^ — M2^- We now construct y in a judicious way in order 
to satisfy the constraint. 

Let y = (0, y^, . . . , Yn) be formed by setting the first vector Y^ equal to 
zero, and, in addition to 1^- (j = 4, . . . , — 2) all lying in LI, we require 
that they are mutually orthogonal of length and are orthogonal to both 
/3 and 7r(Q). This is possible since N < n and can be achieved by Gram- 
Schmidt orthonormalization. Then by construction, Y^Z = () since {Yj,a) = 
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(y„ /3) = (y^, 1) = for J = 3, . . . , iV, and 

/ 1 \ 





Y'Y = plN-2 - 9 



1 



. 



But this is exactly what is required to satisfy the constraint ()40p when 
1(4 = • • • = uat = 0. This solves the lifting problem in the case when 
the rows of A are linearly independent and ^^zg / 0. Note that provided 
7r(a) 7^ 0, the solution given by Y is unique up to orthogonal transformation 
of the subspace of 11 orthogonal to 7r(a) and /3. If 7r(a) = there is a further 
freedom in the construction of Y . 

If now the rows of A are linearly independent and zn = 0, then we 
proceed as above, replacing the condition ^^ai = 1, by ^^ai = 0. This 
simply means that the vector a automatically lies in the plane 11 and there 
is no need to project via vr. Note that the conditions f3£ = J2e o^el^e = 
are still satisfied. Then the construction of the configured star is identical. 

Suppose now that the rows of A are dependent, equivalently, the matrix 
AA^ is no longer invertible. From (|38|) . this is equivalent to ui^ + = 1, 



that is, to U3 = U4 = ■ ■ ■ = UN 
equivalent to the condition: 



n 



0. Then (EOl) and (EU) show that this is 



0. 



j,|2 + (7-2)|^z,| 

e I 

Solutions {zi, . . . , Zn} to ()27p which arise in this case are made explicit in 
Lemma 16.31 

We are now required to solve the system AX = 0, where is the 3x (A^— 2) 
-matrix having zero in each of its entries, subject to the constraint X^X = 
pIn-2- By Lemma ESI we can choose a normalization for {zi, . . . , Zn} such 
that Yle cte = and /?£ = /? is a non-zero constant. We therefore have the 
system: 



ai ■ 


• an 


/? • 


■ p 


1 • 


■ 1 



2^13 



\ 



0. 



\XnZ ■■■ XnN J 

where now Ylii^i = 0. In particular, we require YliXij = for each j = 
3, . . . , N — 2 and the constraint requires that the vectors Xj be orthogonal 
to a, to each other and have length y^. We now proceed in a way similar 
to the non-degenerate case. 

Consider with coordinates (ii, . . . , tn)- Let 11 be the (n— l)-dimensional 
hyperplane ii + • • • + 1„ = as above. Then we are required to find N — 2 
vectors of length in 11 orthogonal to a. The dimension of the space 
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orthogonal to a in 11 is n — 2, so since by hypothesis N < n, this can be 
achieved by Gram- Schmidt orthonormahzation. 

The configured stars obtained by the above constructions are fuh provided 
the configuration matrix W has maximal rank N. This is the case if and 
only if the matrix WW^ given in (j28p is invertible. But a simple inductive 
argument shows that 



\wwH 



p + aui 
(JU2U1 

au]\fUi 



auiU2 
p + au2' 



aUNU2 



aU2UN 
p + crUAT^ 



Recall that since 7 < 1, by Lemma 13^ we have p > 0. Since also p + a > 
(equation ([35]) ). it follows that IPFM/^*] > and WW^ is indeed invertible. 
This completes the proof of the theorem. □ 



Example 6.4. We consider the graph on five vertices below, with solutions 
99 to ([1]) normalized so as to take the value at the central vertex and 1 
on one of the other vertices. The symmetry of the figure means that this 
determines the most general non-constant solution. 



z y 




1 X 



There are two solutions to ([T]) with 7 constant, namely: 

7 = 1/3; X = ±i, y = -1 , z = =Fi; 

7 = 1; x = yz, y = ^±^i, z = 3 ±2V2. 

We reject the latter solution, since the inequality 7 < 1 is violated. Con- 
sider the solution with 7 = 1/3. Let us construct the lift at the bottom 
left-hand vertex. First, we normalize so that the solution takes on the value 
at this vertex: 



—i —1 — i 




1 i + i 
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The choice n = = 3 is determined and from (j34|) and (jSSp . we obtain 
p = 2 and a = 1. From ()30p we find that ui^ + = 1 so that U3 = 0. 
Then the 3x1 -matrix B vanishes and system (|36p has general solution 
X3 = (2A,— A,— A). The constraint ([57|) requires that (X^jX^) = p = 2, so 
that A = ibl/\/3- The star matrix W (whose columns give the positions of 
the external star vertices) is given by: 



W 



( 



\ 



-1 




-1 
-1 

- 1 



V3 "^v^ 

We can proceed similarly with the central vertex of degree 4. Now we can 
choose = 3 or = 4. In either case, lii = M2 = 0, p = 2 and fi = 1, so 
that, for = 3 we must have = ±1. Then 





(l 


-1 0^ 






A = 


1 


-1 








u 1 






^2^3/ 


and the unique (minimizing 


) solution is given by 





A+B 



I 4\ 

2 



The star matrix W is given by 

/ 1 

W = 

\ 2 





1 

V3 



-1 



2 



\ 

-1 

/ 

2 / 



where the last line is only defined up to sign. 

At the same vertex, we can also take A^ = 4. Then 

/ 
B= 

\ 2^/3n3 2^3 

We proceed as in the proof of Theorem 16.11 Take ua = 1 and M4 
solution is then given by the (2 x 4)~matrix X = Z + Y , where 



0. The 



Z = A+B 



I V3 
' 2 

^^ 

2 

vi 

2 

■ %/3 
\ 2 



and Y = (0, Y^, 
ti+t2 + h + t4 = 



, with I4 G a vector of length = \/2 in the plane 
orthogonal to both a = (1, 0, -1, 0) and P = (0, 1, 0, -1). 
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Up to sign, this is given by Y4, = (1/^/2,-1/^/2,1/^/2,-1/^/2). We then 
obtain four possible hfts corresponding to the different choices of square 
root, with star matrix: 



W 



( 1 



2 

V 75 



V2 V2 



-1 

2 
1 

V2 / 



What additional information is required in order to make a unique choice 
for the lifted star? This question is important when we come to consider edge 
curvature in Section [71 For example, given a solution to ([1]) corresponding 
to one of the regular polyhedra, we would like the lifted star to be exactly 
the one that arises from its canonical embedding in Euclidean space. One 
way to do this is to define a notion of orientation as follows. 

In [3] , a notion of orientation was considered on a regular graph of degree 
n, say, whereby the graph is endowed with an edge colouring of the n colours 
{1, 2, . . . , n}. Thus each edge is coloured in such a way that no two edges 
of the same colour are incident at a vertex. This enables one to uniquely 
label the edges at each vertex to give an ordering. One could then attempt 
to apply a right-hand rule say, in order to make a choice of lift. However, 
although this can be done with the solution corresponding to the framework 
of the tetrahedron in a way consistent with its embedding, it turns out to 
be impossible for the cube and the dodecahedron. In the latter examples, 
any edge colouring with three colours leads to at least one of the two choices 
of lifted stars to be directed in the opposite way required. We therefore 
proceed to define an orientation in terms of an edge colouring together with 
an n-form at each vertex of degree n. 

Definition 6.5. Let T = {V^E) he a graph with largest vertex degree equal 
to M . Then an edge colouring of T is an association of one of the colours 
{1,2,... , M,M +1} to each edge so that no two same colours are incident at 
any vertex. By a theorem of Vizing, any graph can he coloured with either M 
or M-\-l colours (see |14] ). We make the convention to choose the minimum 
M colours when possible. Given an edge colouring of T , at each vertex, a 
volume form is an alternating mapping of the edges which takes on the 
value +1 or —1. Thus if x is a vertex with n incident edges ei,...,en 
arranged so that colour fe^^ < colour /e/c J for j < k, then 6x{ei, . . . ,en) = 
±1 with 6'2:(eo-(i), . . . ea-(n)) = sign {a)6{ei, . . . , e^), for any permutation a of 
{1, . . . ,n}. An orientation o/F is given hy an edge colouring together with 
a volume form at each vertex. 
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In order to apply this notion of orientation to make a choice of hfted star, 
in the notation of Theorem 16.11 and its proof, we suppose that = 3 and 
the matrix A of equation ()36p is of maximal rank 3. There is now either 
a unique lifted configured star in the case when U3 = 0, or two choices if 
U3 7^ depending on the sign chosen for = ±^/u^^~+Tl^. The star matrix 
is now given by 



W 



\ Xi3 X23 



\ 



On 
/3n 

Xn3 / 



where the last row is only defined up to sign. Suppose that the vertex 
in question has an orientation according to Definition 16.51 Without loss 
of generality, we can suppose that the edges are coloured with the colours 
{1, 2, . . . , n}, in such a way that the external vertex Xi is joined to the inter- 
nal vertex by the edge with colour 1. Suppose that the volume form satisfies 
9{ei, . . . , e„) = e, where e G {+1, — !}• Then provided the determinant of 
the 3 X 3-minor given by the first three columns of W is non-zero, we choose 
the sign of the third row so that 



Cti OL2 OL3 

Pi P2 /33 

3^13 2^23 3^33 



e5. 



where 5 > 0. If on the other hand this determinant vanishes, then we 
proceed in a lexicographic ordering, to choose next the minor formed from 
columns 1, 2 and 4 and so on, until we encounter a non-zero determinant 
and apply the above rule. 

Example 6.6. Consider the framework of a regular octagon with vertices 
placed at the points (±1, 0, 0), (0, ±1, 0), (0, 0, ±1). Then this can be edge- 
coloured as indicated. Then there is a volume form which gives the lifts 
that correspond to the standard embedding in M^. However, in order to do 
this at the lateral vertices, we have to impose an additional condition that 
the star be regular. This is because at these vertices its = and we do not 
satisfy the conditions of the discussion above. The corresponding solution 
to ([I]) has 7 = 1/2 and p = cr = 2. 
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Consider the vertex xq and define the volume form O^g by 0^q(1234) = — 1 
(for convenience, we write 0(1234) rather than 0(ei, 62, 63, 64)). Then with 
this edge-colouring, at this vertex zi = 1, 22 = i, 23 = — 1 and Z4 = — i. Thus 
ui = U2 = and = ±1. The solution to ()36p is given by = ±(1, 1, 1, 1). 
In order to be consistent with the orientation, we must take the negative 
sign, to give the lifted star: 



W 



\ 



1 


-1 



-1 



-1 



whose sign of the determinant of the 3 x 3-minor given by the first three 
columns is negative, which coincides with the sign of 0a;g(1234). 

At the vertex xi, we choose 0a;i(1234) = +1. Then the edge-colouring 
dictates that = — 1 + i, Z2 = 1 + i, 2:3 = i, Z4 = i, so that ui = 0, U2 = 1 
and ^3 = 0. The solution to (j36p gives a 1-parameter family of lifted stars: 

/ 1 -1 \ 



w 



1 



COS t 

' V2 



cost 
' V2 



-1 





COS t 

V2 



+ sin t 



COS t 

V2 



sint 



If we now impose the condition that the lift must be a regular star, then 
there are just two solutions given by t = 7r/2 or t = 3tt/2. The choice 
t = 3tt/2 is required in order that the determinant of the 3 x 3-minor 
consisting of the first three columns be positive, to coincide with the sign of 
(1234). This gives the lift that coincides with the canonical embedding of 
the octahedron. We proceed similarly with the other vertices, defining the 
appropriate volume form, with the proviso that the the stars at the lateral 
vertices be regular. 



Distance: The above analysis enables us to define edge-length and so dis- 
tance on a graph F = (y,E) admitting a solution to ([TJ, provided that at 
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each vertex we have 7 < 1. For this, note that if we reverse the order of 
multiphcation of W and W^, we obtain 

/ ll^ilP (^1,^2) 



WW 



P2| 



{X2i Xn) 



\ 



\ {xi,Xn) {x2,Xn) 



J 



where {xj^x^) denotes the standard Euchdean inner product of Xi and Xj. 



But then 



El 



trace W^W = trace WW^ 



Np + a\\u\\^ 



Np + a . 



Now this latter quantity can be expressed in terms of 7 and Z£ from ([3; 
and the relation 7 = (t/((T + p) to give the mean of the values \\x(\\: 



(42) 



1 



n 



\Xl)\ 



[N + H- N)-i) 



n(l-7) 

This equation expresses the mean length of the edges of a virtual configured 
star in whose external vertices xi project to z^. This motivates our 
definition of edge length in a graph. 

Let r = {V, E) be a graph coupled to a solution (/? : — C to equation 
([T]). For each x , set 

7(3;) 



p{x) 



(y^x ^ ^ y^x 



where n(x) is the degree of F at x. 

Definition 6.7. If x £ V is a vertex of degree n{x) such that j{x) < 1, 
then we define the median edge length at x relative to 93 to be the quantity 
r{x) whose square is given by 

^2 _ [N + {I - N)-f{x^^ 



r{xy 



p{x) . 



n{x)[l — 7(x)] 

Ifxy £ E is an edge which joins x to y such that both 7(3;) < 1 and 7(7/) < 1, 
then we define the length of xy relative to ip to be the mean i{xy) of the 
median edge lengths at x and y: 

iixy) = !M±lM . 

As emphasized in the above definition, the lengths so defined are relative 
to the solution of ([1]), which is only defined up to 99 1— t- \ip + p for A, G C. 
This means that the only meaningful quantities are relative lengths, say 
£{e)/£{f), for two edges e, / G E. This is consistent with our relational 
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interpretation of physical quantities as discussed in the Introduction. In 
particular, if both n and 7 are constant on the graph, we may take the 
quantity 2p defined by as a measure of median edge length at each 
vertex: 

r{xf = 2p = Y, W{y) - ^{^)? - 1\ Y.^^{y) - v^(x))f . 

y^x y~x 

We can define an absolute length by normalizing as follows. Let V = 
{V^E) be a graph coupled to a solution 99 : 1/ — )• C to ([1]). Then as in 
Appendix [Cl define the square L'^-norm of the derivative of ip to be the 
quantity: 

xy&E x,y(^V,x^y 

where df^xy) = ip{y) — ^{x) is the discrete derivative with respect to some 
orientation of the edge xy (in this case x is the initial vertex and y the end 
vertex) . 

Definition 6.8. Let T = {V, E) he a graph coupled to a solution 99 : 1/ — t- C 
to equation ([T|). If x ^ V is a vertex such that j{x) < 1, then we define 
the absolute median edge length at x relative to ip to be the quantity Tabs (2;) 
whose square is given by 

2 _ ^(a^)^ 

where r(x) is the median edge length at x relative to p. If e £ E is an 
edge joining x to y such that both j{x) < 1 and 7(7/) < 1, then we define 
the absolute length of e relative to ip to be the mean ^abs(e) of the absolute 
median edge lengths at x and y: 

„ , X ?'abs(a;) +?'abs(y) 

4bs(e) = ^ . 

Then both the quantities r^^{x) and ^bsl^j) are independent of the free- 
dom (/? I— + /i (A, /X G C). 

The median edge length of Definition 16.71 is defined so as to give the 
length of the edges of a corresponding regular star in M^, when such ex- 
ists. In particular, if F = (1/, E) is the 1-skeleton of a regular polytope and 
(/? : y — )• C associates to each vertex its value after an orthogonal projection, 
then the edge-length at each vertex coincides with the lengths of the edges 
of the regular polytope. More generally, we can interpret the edge length at 
each vertex as the length of the edges of the "best fit" polytope at that ver- 
tex. The median edge length then gives the average length at two adjacent 
vertices. 
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Example 6.9. If we return to Example 16.41 and consider the solution cor- 
responding to 7 = 1/3, then the (unique) common dimension to define edge 
length is = 3. At the central vertex the edge length is \fl /2 and at 
any of the other vertices, it is y^7/3. Thus the median edge length of the 



edge joining the central vertex to one of the other vertices is ^ ' 

whereas the median edge length of one of the outside edges is y^7/3. So, 
for example, the shortest path joing a; to z is given by passing through the 
central vertex. Note that, as already remarked, the edge lengths are only 
defined up to multiple and so only relative edge lengths have meaning. 

A question we now consider, is whether the notion of distance, either 
relative or absolute, that we have defined above, endows a graph with the 
structure of a •path metric space in the sense of M. Gromov [18j. We first of 
all note a triangle inequality around complete subgraphs on three vertices. 

Given a function ip : V ^ C and a vertex x €V, we say that (p is constant 
on the star centred on x if the restriction of to x and its neighbours y ~ x, 
is constant. 

Proposition 6.10. (Local triangle inequality) Let T = {V,E) be a graph 
coupled to a solution ip : V ^ C to equation ([T|). Suppose x,y,z G V are 
three vertices of a complete subgraph: x y, y z, z x, such that the 
inequality j < 1 is satisfied at each vertex. Then the triangle inequality is 
satisfied: 

£{xy) + £{xz) > l{yz) . 

If further ip is non-constant on the star centred on x, then the inequality is 
strict. 

Proof. This is an immediate consequence of the definition. Specifically, 

i{xy) + tixz) = ^(r(x) + r{y)) + ^(r(x) + r{z)) = £{yz) + r{x) > £{yz) , 

since because of the inequality 7(x) < 1, we have r{x) > 0. If further, ip is 
non-constant on the star centred on x, then r(x) > and the inequality is 
strict. □ 

In spite of this local triangle inequality, we may encounter a difficulty 
in trying to endow a graph coupled to a solution (/? to ([1]) with a metric 
space structure. This may arise when, for a given vertex x, the function ip is 
constant on the star centered on x, as well as on the star centred on one of its 
neighbours y. Then £{'xy) = 0. We can either agree to allow distinct points 
to have zero distance between them, and so consider rather a pseudo-metric 
space structure, or we can avoid this situation by introducting a notion of 
collapsing. This is a concept we will return to in Section [9l 
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Definition 6.11. Let {T,ip) be a graph coupled to a solution to equation 
([1]). Then we collapse T to a new graph T by removing all edges that con- 
nect vertices at which f takes on identical values; then remove all isolated 
vertices. 

It is clear that after collapse, if we let if denote the restriction of to 
r, then if also satisfies ([TJ with 7 = nj/n where n is the new degree at 
each vertex. Indeed, if we check at a vertex x, then if y is a neighbour at 
which ip{y) = ^{x), then since only the difference ^p{y) — <p{x) occurs in ([1]), 
removing the edge xy only affects the degree. However, it is to be noted 
that collapsing may disconnect a graph. 

Let {T,(p) be a graph coupled to a solution to equation ([TJ. Then given 
a path c := X0X1X2 • • • Xp joining x to y (so we have x = xq, y = Xp and 
Xj ~ Xj^i for all J = 0, ... ,p — 1), then we define the length i(c) to be the 
sum: 

p-i 

^£{xjXj+i). 

j=0 

We can now define the distance between two vertices to be the infimum of 
the lengths of all paths joining the two vertices. Then provided T is collapsed 
with respect to 93, it is clear that this notion of distance endows T with the 
structure of a path metric space. 

7. Curvature 

Our introduction of curvature on a graph is based on Theorem 16.11 and 
Corollarv 14.21 Thus, we consider a graph T = {V, E) coupled to a solution Lp 
to equation ([1]): 7(^(A(/?)^ = dip^ . At each vertex y G we measure the (in 
general solid) angular deficit 5{y) as determined by a regular star placed in 
whose central vertex xq has the same degree as y in F, and which also 
solves ^ at xq with 7star(2;o) = 7</j(y)- 

Depending on the degree n and the value of 7, 5{y) will be well-defined. 
However, in some situations, there may be different possibilities for the 
dimension N ^ leading to different possible values for the curvature. The 
restriction 7 < 1, will be a necessary condition. 

The aim is to suppose 7 is part of the geometric spectrum (see Section 
[8]), so that both dimension and curvature arise from purely combinatorial 
properties of the graph (independent of (p). For some graphs, these will be 
uniquely defined. This is the case for the bipartite graph -fCas, for example, 
whose geometric spectrum contains the unique value 7 = 1 (see Section [8]) 
and whose vertex degree dictates that it "lives in " dimension three. 

In what follows, we discuss convex polytopes, which are by definition, the 
closed intersection of half-spaces (whether this be in Euclidean space, or in 



50 



PAUL BAIRD 



spherical space). In the case when a polytope is regular (convex or not), its 
vertices all lie on a sphere called the circum- sphere [13]. It is useful to use 
absolute angle measure when measuring solid angles (see [291 |T9]). We will 
write H^'^ (A) for the M-dimensional Hausdorff measure of a set A in these 
units. Then, in any dimension, the angle is measured as a fraction of the 
total angle subtended by a sphere centred at the point in question. Thus in 
two dimensions, a right-angle has value 1/4, whereas in three dimensions, the 
angle subtended by the vertex figure of a cube has value 1/8. Equivalently, 
H'^{A) = 1/8, where A C S'^ represents one eighth portion: x,y,z > 0, of 
the sphere + + = 1 . 

Definition 7.1. (A^-dimensional vertex-curvature) Let{T,(p) be a pair con- 
sisting of a graph T coupled to a solution if to Let y be a vertex of T and 
let n be the degree of T at y. Then the N -dimensional vertex- curvature at 
y is defined provided 'y{y) < 1 and there is a regular N -polytope with vertex 
figure P (a regular (N — l)-polytope) having n vertices. Let P be centred on 
G M^"-*^ with vertices vi, . . . ,Vn lying on its circum- sphere of radius r > 0. 
For 7 < 1, let U be the corresponding configuration matrix with associated 
constant p> (see ([7|) and Let 



be the corresponding vertices of a regular star in centred on 0. Let A 



in absolute angle measure. In the case when 7 = 1, then we define the N- 
dimensional vertex-curvature to be the limit lim^_^i- S^{y), when this exists. 

Note that in the above definition, if we rescale the vertices V£ by Vi i— )■ 



curvature S^{y) are well-defined and independent of this scaling. We write 
vertex-curvature to distinguish it from edge- curvature which we define later, 
but if the context is clear, we shall just write A^-curvature. 

We first of all justify this definition and put it into the context of classical 
work on topological invariants of polytopes. In particular, we refine the 
definition for 3- and 4-curvature. We begin by reviewing a generalization 
of a theorem of Descartes, by Shephard [29], Griinbaum and Shephard [19] 
and Ehrensborg [16j. The context is that of elementary polytopes. 

A family {Fi, . . . , Fr} of (A^— l)-dimensional convex polytopes in form 
an elementary polytope P of dimension N if: (i) for all j, k, Fj D F^ is either 





say, then r i— t- r' 



Ar, p p' 



}?p and both Xj and the 
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empty or a face of each of Fj and F/.; (ii) Uj-Fj is an (A^ — l)-dimensional 
manifold. 

Given an elementary polytope P C M^, denote by V the face decom- 
position of P; thus V is the collection of all (open) faces of all dimen- 
sion, consisting of the vertices, edges, ... , [N — l)-faces, A^-faces. For 
X € M^, w S S'^~^, following Ehrensborg [16], we define the quantity 
R{x,w) := lim5^o+ ^pi^ + s ■ w), where Ip is the characteristic function 
of P. Note that this takes on the value or 1. Then given a face F € V, we 
have R{x, w) = R{y, w) for all x,y £ F; write R{F, w) for this and define 
AF = {we S^-^ : R{F, w) = 1}. 

Let S{Af) = H^-^idAp) and let aN = {S^), so that in absolute 
angle measure, = 1. Let F be a face of P of dimension < — 3. We 
define the deficiency at F to be the quantity: 

6{F) := aN-2 - S{Af) . 

Note that if Q C S^~^ is spherically convex (that is, it is the intersection of 
hemispheres), then S{Q) = H^~'^[dQ) is proportional to the Haar measure 
of all the great circles which intersect Q. The following theorem generalizes 
a classical result of Descartes. 

Theorem 7.2. [29^ [T9l [T6] Let P he an elementary polytope with face de- 
composition V such that P has only one N- dimensional face P^ . Then 

e{F)6{F) = a^^,[{-lf - l]e{P') , 

FeV, dim F<N~3 

where £{F) denotes the Euler characteristic of F given by (— l)'^ when F is 
of dimension k. 

In the case when = 3 and P is a convex polyhedron (now using radians 
for our measure), we obtain the classical theorem of Descartes: 

'^Sivi) = Att, 

I 

where the sum is taken over the vertices of P . Here, the deficiency 5(u£) is 
the sum of the internal angles at vi of the faces which contain vi. We may 
view this theorem as a discrete version of the Gauss-Bonnet Theorem for 
surfaces in the smooth setting. 

Let us now consider the different dimensional curvatures that arise from 
Definition 17.11 By convention, at a vertex of degree 1, we assign the cur- 
vature (5 = 1. It is straightforward to see that at a vertex of degree 2, the 
2-curvature just measures the exterior angle in absolute angle measure. In 
view of identity (l6|), we can state this as follows. 
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Proposition 7.3. Let (F, ip) be a pair consisting of a graph T coupled to 
a solution ip to Let y be a vertex of degree 2 where 7 < 1, such 

exists. Then the 2 -dimensional vertex-curvature of (T, (p) at y is given by 
the quantity: 



Note that lim^_^x- ^-yiu) is well-defined and equals 1/2. For example, if 
r is a cyclic graph of even order 2k and (/? is a function taking on alternate 
values at neighbouring vertices. Then ip satisfies ([1]) with 7 = 1. The total 
curvature is then given by A;. If T is a regular polygon in the plane and (/? 
the corresponding position function, then the total 2-curvature is equal to 
1, or in radians, to 2tt, as required. 

We now proceed to higher dimensional curvature; dimension 3 is of par- 
ticular interest because of the minimizing property of the solution to the 
lifting problem that occurs in this case; see Theorem 16.11 

Proposition 7.4. (3-dimensional vertex-curvature) Let (T, ip) be a pair con- 
sisting of a graph T coupled to a solution ip to Let y be a vertex of T and 
let n be the degree ofT at y. Suppose that n G {3, 4, 5} and that 7 < 1. Then 
the ?>- dimensional vertex- curvature of (F, ip) at y is given by the quantity: 



(£ = 1, . . . , n). The boundary of the convex hull A of the set {xi, . . . , x„} in 
5^ is made up of arcs of great circles of length a. = arccos(xi-X2)- In absolute 
angle measure, the deficit, or 3-curvature, is given by 1 — Substitution of 
the expressions for xi given by Definition 17.11 gives the required formula. □ 

Our requirement that n G {3, 4, 5} is a consequence of Definition 17. H for 
the only polygons that appear as vertex figures of regular polyhedra have 
these possibilities for their number of sides. Of course, the expression for 
6^{y) above is defined for any n provided 7 < 1. 

Note that \\m^^i- 6^{y) is well-defined and equals 1. For example, the 
bipartite graph ^"33 has the unique value 7 = 1 in its geometric spectrum. 
Also the degree of each vertex is n = 3. The 3-curvature at each vertex 
is then 5 = 1 and the total curvature is given by X^ygy'5(y) = 6. On the 
other hand, the 1-skeleton of a tetrahedron in has n = 3 and 7 = 3/4, so 
that 7 < 1. It is easily checked that the 3-curvature at each vertex is given 
by 5 = 1/2, giving a total curvature of 2. To obtain the total curvature in 
radian measure, we multiply by 27r to give the value 47r, which confirms the 
theorem of Descartes. 





Proof. The configuration of vertices V£ is given by ()13p . that is ve = e 



.2nie/n 
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For degree 3, the 3-curvature is the only A^-curvature that can apply, 
since there is no other vertex figure with three vertices. This then gives a 
well-defined curvature for vertices of degree 3, provided 7 < 1. 

Example 7.5. Consider the double cone on the triangle discussed in Section 
[5j The corresponding framework satisfies equation ([1]) with 7 = 4/5 on the 
three lateral vertices of degree 4 and 7 = 1/3 at the two apexes of degree 3. 
Then we calculate: 

Oanex = 1 arccos - and Oiat = 1 arccos - , 

apex 7 2tt 7 

to give total curvature: 

3/1 5 
<5tot = 3(5iat + 2(5apex = 5 - - ( arccos J arccos - 

Now 



cos 



^arccos ^ + 2 arccos = ^{l - 240\/2) = -0 • 98662 . 
to five decimal places, so that 

arccos — \- 2 arccos — ~ vr . 
7 7 

is close, but not equal to vr. The exact value vr gives a total curvature of 
2, which is the value that we expect from the theorem of Descartes. The 
small difference arises due to the fact that the vertex figures at the lateral 
vertices are not configured stars, whereas the curvature is defined in terms 
of the deficit that occurs for the unique lifted configured star. That is, we 
try to fit a regular polytope in the best way possible. What is remarkable, 
is how close the two values are. 

The above example illustrates one of the problems in defining the cur- 
vature. The advantage of lifting to a configured star is that, in dimension 
N = 3, the lift is unique and so the curvature is uniquely defined. However, 
any expression of the total curvature as an invariant quantity would need to 
involve some approximation. 

For the 4-curvature, there are some special cases to consider. We list 
these in the proposition below. 

Proposition 7.6. (4-dimensional vertex-curvature) Let (P, ip) be a pair con- 
sisting of a graph T coupled to a solution ip to Let y be a vertex ofT 
and let n be the degree of T at y. Suppose that n G {4,6,12,20} and that 
7 < 1. Then depending on the degree, the 4- dimensional vertex- curvature of 
(r, if) at y is given by one of the expressions below: 
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degree vertex figure 



curvature 



12 



20 



tetrahedron 



octahedron 



icosahedron 



dodecahedron 



3 / 7 
arccos 

6 / 1 ' 

arccos 

vr \ 5 — 37 ^ 

/ 6(\/5 + 1)7 - 11 - 7V5 

— arccos ^ ^ 

TT \ 2[6(\/5 + 3)7- 23- 7\/5] 



15 



57 - 1 - 2V5 



15 

10 arccos , ^ 

vr \ 2[-57 + 8- ^/5] 



Proof. Given two vectors u,v ^ 5*^(r) in a sphere of radius r, the arc of the 
great circle joining n to iT is given by 



cos ( 



[u ■ V) sin( 



\Jr^ — (u ■ vy 



u + 



sin ( 



(0 < < arcsin 



■\J r'^ — (u-v)'^ 



When r = 1, this is unit speed. Furthermore, the tangent vector to this arc 
at V is given by 

1 



(43) 



===(— r M + (n • v)v) . 



The area of a spherical polygon with m sides and with interior angles 6k 
(fc = 1, . . . , m) is given by 

KA = ^ek- {m-2)'K , 



where K is the curvature of the sphere. We are required to calculate the 
spherical surface area of the boundaries of the various vertex figures in S'^ . 
These are made up of faces lying in great 2-spheres which are either triangles, 
or in the case of the dodecahedron, pentagons, whose edges are arcs of great 
circles. In order to calculate the interior angles, we calculate the scalar 
product of the unit tangents to these edges at a vertex. By symmetry, any 
vertex will do. We calculate this for the icosahedron and the dodecahedron, 
the other cases being similar. 

For the dodecahedron, a configuration of vertices is given by (jlSp . With 
the notation of Definition 17. H three consecutive vertices around one pentag- 
onal face are given by 

/ a \ 



1 



Xl 



V8 + 3a2 



a\ 


V V8 / 



1 



X2 



V8 + 3a2 



a 
a 
a 

\V8j 



1 



X3 



V8 + 3a2 



\ 



/ 

aX 
\ V8 J 
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where a = ■\/n{l — 7) and A = (1 + \/5)/2. These three vertices determine a 
great 2-sphere in which contains the pentagonal face. With this arrange- 
ment, X2 is the central vertex joined to xi and X3 by arcs of great circles. 
In order to calculate the interior angle at each vertex of the pentagon, we 
calculate the tangent to each of these arcs at X2- To do this we apply (jl3|) . 
For the first arc we set u = xi and v = X2, to obtain the tangent vector: 

/ 4(3-^/5) + 4(3-^/5) \ 



1 



2V3a2 + 8\/a2 + 12 - 4^5 



-^(3 + ^/5)+4(l-^/5) 
Vda"^ + 8 



\ 



For the second arc, we set u 



t2 



2V3a2 + 8\/a2 + 12 - 4\/5 



-a\/8(3 - Vb) 

X3 and V = X2, to obtain: 

/ ^/5a2 +8 \ 

4(3-^/5) + 4(3-^/5) 
-4(3 + ^/5) + 4(1-^/5) 



Then 



tl-t2 



\ -aV8{3 - Vb) 

-a? + 16 - 8^/5 



/ 



2(o2 + 12 - 4^/5) ' 

which gives the cosine of the interior angle (it is indeed the interior angle, 
being greater than vr/2). Then the area (in absolute angle measure) of each 
pentagonal face is given by 

-0^ + 16-8^/5' 



5 arccos 



37r 



2(a2 + 12-4V5)^ 

so that the surface area of the spherical dodecahedron is given by twelve 
times this quantity. We then obtain the angular deficiency, or 4-curvature: 



10 



15 



arccos 



-a2 + 16 - 8\/5 \ 



TT 



2(a2 + 12 - 4^/5) j 

On substituting the value of a, we obtain the required formula. 

For the icosahedron, a configuration of vertices is given by (jl7p . Three 
vertices which form one of the triangular faces are given by: 













( ^ 


Vi = 


(') 


, V2 = 


(') 













^0/ 




^ 1 



where A = (1 + \/5)/2. Then p = 2 + = b + y/b, n 
(5 + \/5)/2. This gives the corresponding vertices in S as: 



12 and r 



xe 




^/2 + n{l 



1,2,3). 
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We then proceed as above for the dodecahedron to calculate the angular 
deficiency. □ 

Example 7.7. The 600-cell is a convex 4-dimensional regular polytope made 
up of 600 tetrahedral 3-polytopes. It has 120 vertices and 720 edges. Its 
vertex figure is a regular icosahedron. If we consider an orthogonal projec- 
tion onto the complex plane and let (p associate the corresponding value to 
each vertex, then by Theorem ip satisfies ([1]) with 7 constant. We can 
find the value of 7 as follows. 

Since the edges of the 600-cell all have the same length, in the notation of 
the above proof, we must have the distance from the origin to xi, that is 1, 
equal to the distance between two neighbours of the vertex figure: | |xi — X2 1 1 . 
One can readily calculate: 

11^ ^1,2 8^(1 - 7) 

KCi — xo = ^ , 

(5 + V5)[2 + n(l-7)] 

to obtain the negative value: 

5(1 - 2^/5) 
^ = 3 • 

One can now confirm the generalization of the theorem of Descartes (The- 
orem [721) • 

The 600-cell has 5 tetrahedra around each edges, each having dihedral 
angle arccos(l/3). Thus the angular deficiency at each edge (in absolute 
angle measure) is given by: 

Op = 1 arccos — . 

27r 3 

Substitution of the above value of 7 into the third formula of Proposition 
17.61 gives the deficit, or curvature at each vertex, as 

0^ = 6 arccos - . 

vr 3 

On then finds that 120(5^, — 720 = 0, as required. 

We can proceed similarly to obtain explicit formulae for higher dimen- 
sional A^-curvature. This is simplified by the fact that there are just three 
regular polytopes in dimensions > 5, namely the A^-simplex, the A^-cube 
and the cross-polytope, with vertex figures an (A^— l)-simplex in the first two 
cases and another cross-polytope in the last case. To find the A/'-curvature 
requires the calculation of the (A^ — 2)-dimensional measure of (A'^ — 2)- 
simplices in great spheres of S^~^. This is a standard, but non-trivial 
procedure using Schlafli's differential equality [28]. See also the expository 
article of J. Milnor for a nice account and references pi]. The article of 
J. Murakami provides explicit expressions in the 3-sphere j22j . We do not 
attempt to derive these formulae here. 
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Edge- curvature Let F = (V, E) be a graph endowed with a solution (/? to ([1]), 
together with a choice of hft of a configured star into at each vertex, 
where the dimension is to be fixed over the whole graph. If we suppose 
= 3, this may be achieved by defining an orientation on F, as discussed in 
Section [6j We suppose further that each star has a well-defined axis defined 
by a unit vector u{x) S M^, for each x £ V . The axis should be directed 
from the internal vertex of the star towards its centre of mass. 

Definition 7.8. Given an edge e = xy G E, define the edge-curvature of 
e to he the unique angle 9{e) := arccos{{u{x),u{y))^N) G [0,7r]. Given a 
vertex x £ V , define the mean edge-curvature at x to be the mean of the 
edge-curvatures of the edges incident with x. 

Thus the edge-curvature measures the angle between the axes of adjacent 
stars. It is clearly independent of the freedom i— t- XLp-\- fj, in the solution (p. 
By analogy with Riemannian geometry, various other curvatures can now 
be defined. If we let £{e) denote the length of an edge e = xy as given by 
Definition 16.71 and 0{e) its edge-curvature, then the radius of the best-fit 
circle is given by r(e) = l{e)/9{e) (by best-fit circle, we mean the circle 
subtending the same arc length i{e) for the given angle 9{e)). The normal 
curvature of e is then the reciprocal l/r[e) = 9{e)/i{e). The mean curvature 
at a vertex x is the mean of the normal curvatures of the edges incident with 
X. Since i{e) depends on the scaling ip i— )• Xip, this quantity also depends on 
the scaling; the mean curvature should be thought of as the analogue of the 
same notion in the smooth setting, wherby we locally embed a Riemannian 
manifold in a Euclidean space. 

Ricci curvature is one of the most natural curvatures intrinsic to a Rie- 
mannian manifold. Recall that given two unit directions X and Y, the 
sectional curvature Sec{X ,Y) can be interpreted as the Gaussian curva- 
ture of a small geodesic surface generated by the plane X AY. This in 
turn is the product of the principal curvatures which are the extremal val- 
ues of the normal curvatures. The Ricci curvature Hie {X , X) is then the 
sum: ^jSec {X,Yj) taken over an orthonormal frame {Yj} with each Yj 
orthogonal to X. This motivates the following definition. 

Definition 7.9. Given a vertex x £ V and two edges ei = xyi and 62 = xy2 

with endpoint x, we define the sectional curvature Sec (ei, 62) determined by 
ei and 62 to be the product: 



Secx (61,62) = 9{ei)9{e2) , 
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where 9{ej) (j = 1,2) are the edge- curvatures. For an edge e = xy, the Ricci 
curvature Ric (e, e) is the sum 

Ric,(e,e) =^(e)2 ^ e{xy)e{xz) , 

and the scalar curvature at x is given by 

ScaU = ^Ric {xy,xy)/i{xyf. 

The length scaling does not appear in the sectional curvature, since in 
the smooth setting this quantity depends only on the plane generated by 
two unit vectors. On the other hand, the Ricci curvature is bilinear in its 
arguments and so should depend on the square of the length. In Riemannian 
geometry, one usually applies the polarization identity to define Ric {X,Y), 
however, there would seem to be no reasonable interpretation for the sum of 
two edges in our setting. The dependence on length is once more removed 
from the scalar curvature, which is the trace of the Ricci curvature. 

8. The geometric spectrum and the 7-polynomial 

Recall the geometric spectrum of a graph F = {V, E) is the set 

S = {7 S M : 3 non — constant ip :V ^ C such that 7(A93)^ = (d^?)^} . 

It is clear that the spectrum only depends on the isomorphism class of a 
graph. Section O shows how a particular value in the spectrum may corre- 
spond to local Euclidean geometry. Thus, at a particular vertex, the edges 
which connect it to its neighbours may be realised as vectors in a Euclidean 
space; in particular their relative lengths are defined as well as the angles 
between them. The edge-curvature as defined in Definition 17. 8| may then be 
considered as a measure of how these local Euclidean geometries are pieced 
together to form a global geometric object. None of these aspects require 
the graph to be embedded in an ambient space; they emerge purely from 
the combinatorial properties of the graph. 

There are some obvious questions about the geometric spectrum of a 
graph: is it discrete? is it finite? are there bounds? To compute it, even for 
simple graphs, is quite challenging. We deduced in Section [3] the spectrum 
of some cyclic graphs of low order. However, once the order increases, then 
the problem can become difficult. Real regular cyclic sequences correspond 
to polynomial equations over the integers with positive coefficients having 
real roots; complex solutions correspond to closed walks in the plane with 
each step forming an angle ±0 with the previous step, for some fixed 6. We 
may attempt an algebraic geometric approach to shed some light on these 
issues. 
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Let r = {V, E) be a connected graph. We are interested in the possible 
real numbers 7 for which there are non-constant solutions to the equation: 

Any solution is invariant by (/? 1— t- Xip + /i, for complex constants A, /i. Con- 
sider first how to parametrize all possible complex fields on the graph under 
this invariance. 

Label the vertices of the graph xi,X2, ■ ■ ■ ,X]\f and consider a non-constant 
complex field <p that assigns the value ip{xk) = Zk to vertex x^- Then the 
space of all such fields is identified with the complex space C^\{fi{l, 1, . . . , 1) : 
/i G C}. Up to the equivalence (zi, . . . , zn) ~ (-Zi + /x, . . . , ztv + /^), we can 
identify these fields with the set H \ {0}, where n is the linear subspace 
n = {Z = (zi, . . . , zat) G : zi + • • • + z„ = 0} C C^. In effect, given 
any non-constant field (zi, . . . , z^r), then (zi + /x, . . . , Zjv + /i) lies in the 
plane zi + • • • + zat = 0, when we set /x = — -^(zi + • • • + zj^). By non- 
constancy, this is non-zero. Furthermore, it is clear that any two equivalent 
fields correspond to the same point. 

Now consider the equivalence Z ~ \Z, for A G C \ {0}. This defines 
the moduli space of fields up to equivalence, as Z := CP^'"^. Specifically, 
given a point [zi, . . . ,ziy-i] £ Z in homogeneous coordinates, we define a 
representative field by (zi, . . . , zn-i, zn = — Ylk=i ^k) £ C^. In practice, 
we can set a field equal to and 1 on two selected vertices xq and xi, 
respectively, and label the other vertices arbitrarily. This is only one chart 
and we miss those fields which coincide at these two vertices. 

If we consider 7 as an arbitrary complex parameter, then ([T]) imposes a 
constraint at each vertex, so we have equations in A — 1 parameters. 
In general these are independent so that this is an overdetermined system, 
which may have no solutions. The graphs on five and six vertices below have 
empty geometric spectrum. 




Figure 2. Two graphs which admit no non-trivial solutions to ([T|) with 7 

constant. 

Even with such simple examples, the equations are quite difficult to solve 
by hand. Let us consider one way to approach the problem of computing 
the geometric spectrum. 
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As above, let F = {V, E) be a finite connected graph with N vertices 
labeled xi, . . . , xn- For each I = 2, . . . N, consider the following set of N 
polynomials defined over the algebraically closed field C. The variables 
are the values {zi, . . . , z^} of a field on F with constraints zi = and 
Zi = 1; we suppose the degree of vertex j is n{j) and that zj^ £ {zi, . . . , zn} 
{k = 1, . . . , n{j)) are the values of the field on the neighbours Xjk of xj. The 
polynomials are then defined by 

\ ^ ri(j) 

^{zj - Zjk) - ^{zj - Zjkf {zi =0,ze = l), 
k=l J k=l 

in the — 1 complex variables {7, Z2, z^, . . . , zi, . . . , Z]\f}. Recall some facts 
and terminology from commutative algebra. We are particularly interested 
in the techniques of Grobner bases, for which we refer the reader to [HET]. 

For an ideal / =< /i, . . . , /at > in a polynomial ring C[xi, X2, • • • , xm], 
we denote by V{I) the corresponding variety given as the solution set of the 
equations /i = 0, /2 = 0, . . . , /jv = 0. Then I is called zero- dimensional if 
V{I) is finite. A Grobner basis for / is a basis of polynomials which can be 
constructed from fi, ■ ■ ■ , /n using a particular algorithm, called the Buch- 
berger algorithm. To employ this algorithm, one is required first to choose 
an order on monomials. We shall only be concerned with lexicographical 
order here, which means we first choose an ordering of the variables, say 
xi > X2 > • • • > Xm and then order monomials := xi°^ • • • XAf°", 
:= xi^i • • • xm'^*^) by x" < x^ if and only if the first coordinate ai and 
/3j from the left which are different satisfy Ui < Pi. With respect to the 
monomial order, every polynomial f in I has a leading term It (/) which is 
the product It (/) = Ic (/)lm (/) of the leading coefficient with the leading 
monomial. 

A set of non-zero polynomials G = {gi, . . . ,gp} in / is called a Grobner 
basis for I if and only if for all / G I such that / 7^ 0, there is a gj in G 
such that Im (gj) divides Im (/). The Grobner basis is further called reduced 
if for all j, Ic (gj) = 1 and gj is reduced with respect to G \ {gj}, that is, 
no non-zero term in gj is divisible by any Im (^r^) for any k ^ j. A theorem 
of Buchberger states that every non-zero ideal has a unique reduced Grober 
basis with respect to a monomial order [TJ. Grobner bases are particularly 
useful for understanding the solution set of a system of polynomial equations. 

Let I be an ideal in the polynomial ring C[xi,X2, • • • ,xm] and let G = 
{gi, . . . ,gp} be the unique reduced Grobner basis with respect to the lexi- 
cographical ordering induced by the order xi > X2 > • • • > x^f . Then V{I) 
is finite if and only if for each j = 1, . . . , M, there exists a, g^ £ G such 
that Imgifc = Xj"^ for some natural number Uj. As a consequence, if / is 
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a zero-dimensional ideal, it follows that we can order gi, . . . ,gp so that gi 
contains only the variable xm, 92 contains only xm,xm-i and so on. This is 
because the leading monomial of one element, gi say, of G must be a power 
of X M and then no other term of gi can contain powers of any other variable 
(for such terms would be greater that any power of xm with respect to the 
monomial order), and so on for successive elements g2,g3, ■ ■ ■ of G. We note 
also that V{I) is empty if and only if 1 G G. 

It is also the case that, with the above hypotheses, the polynomial gi is 
the least degree univariate polynomial in xm which belongs to I (any zero- 
dimensional ideal contains such a polynomial for every variable). For if 
there was another univariate polynomial p{xm) with degp < deg^i, then 
Imp would divide Imgi in a strict sense, which would contradict the fact 
that G is a reduced Grobner basis. Let us now return to the case under 
consideration. 

For each £ = 2, . . . , N , consider the ideal Ii =< fi^, . . . , Jn^ >. Suppose 
that for each i = 2, . . . , N this admits a least degree univariate polynomial 
Pi in ^. This can be constructed by first choosing a lexicographical ordering 
of the variables with 7 the smallest and then applying an algorithm (say the 
Buchberger algorithm) to construct the unique reduced Grobner basis for 
The first element of this basis gives pi. 

Definition 8.1. We define the 7-polynomial p = pr of the connected finite 
graph T = (V, E) to be the least common multiple of the least degree univari- 
ate polynomials pi (i = 2, . . . , N ) in ^ associated to the equations ([T|) for 
fields (zi, . . . , zjv) on T with zi = and Z£ = 1: 

p := 1cm (p2, • • • ,Pn) , 

when each pi exists. 

The 7-polynomial ^(7) is defined up to rational multiple and has rational 
coefficients. This is because the initial polynomials fj^ used to define p all 
have integer coefficients and the Buchberger algorithm then generates poly- 
nomials with rational coefficients-it involves at most division by coefficients- 
see [1]. Clearly p depends only on the isomorphism class of a graph and in 
the case when the equations ([T|) admit no solutions for 7 constant and com- 
plex, then p = 1. In this case we shall say that p is trivial. The polynomials 
Pi and so p may still be well-defined even if the solution set of the equations 
is infinite (that is the corresponding ideal is no longer zero-dimensional) . In 
fact we know of no case when they are not well-defined. 

The elements of the geometric spectrum arise as real roots of p (the prob- 
lem of establishing the discreteness of the spectrum is clearly intimately 
related to knowing if p is well-defined in all cases). However, not all real 
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roots may occur in the spectrum, for in general they must also solve the other 
equations determined by the Grobner basis: gi = 0, . . . , gp = 0. Examples 
below illustrate this property. We know of no two non-isomorphic connected 
graphs with non-trivial 7-polynomial having the same 7-polynomial. How- 
ever, the examples of Figure 2, give two non-isomorphic graphs having trivial 
P- 

The examples of the triangle C3 (the cyclic graph on three vertices) and 
the bipartite graphs K23 and are instructive. We label the vertices as 
indicated and consider fields if taking the values ip{xj) = Zj at each vertex 




Xl X2 Xl X2 Xl X2 X3 

C3 K23 K33 



For the triangle, there are precisely two solutions to ([T]) when we normalize 

1 /s 

so that zi = 0,Z2 = 1; specifically z^ = ± i^. Then p = p2 = ps = 3j — 2 
is the 7 polynomial and the geometric spectrum is the unique root 7 = 2/3. 

For K23, we find P2 = 1 with no solution and p3 = 7^ — 27 + 1 = (7 — 1)^ 
with solution zi = 0, 23 = 1, = 0, 25 = A arbitrary and Z4 = [1 + A ± 
\/3(l — A)i]/2. Then p = 7^ — 27 + 1 and the geometric spectrum is given 
by S = {1}. 

For K33, we find p2 = 97^ - 267 + f 7 = (7 - 1)(97 - 17) and p^ = 9-/^ - 
357^+437—17 = (7—1)^2, so that the 7-polynomial p = 97'^— 357^+437— 17. 
Although this has 7 = 17/9 as a root, the geometric spectrum S = {1}. In 
fact for 7 = l,zi = 0,Z2 = 1 we find a two complex parameter family of 
solutions as in Section [5J The next example shows that even for simple 
graphs, the 7-polynomial can be quite complicated. 

Consider the graph of constant degree three on six vertices whose edges 
form two concentric triangles as shown below. 
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The 7-polynomial is given by 

58593757^° - 676562507^ + 3335218757^ - 9260250007^ 
+16039788307*^ - 18084860287^ + 13396555987^ 
-6398928727^ + 1867603237^ - 295988587 + 1883007. 

This has eight real roots, four of which are rational: 7 = 3/5,21/25,1,3. 
The value 3 lies in the spectrum and corresponds to the obvious colouring of 
the vertices with two colours, by choosing identical colours for each triangle. 
The value 1 also lies in the spectrum and corresponds to the colouring of 
each triangle with the three colours 0, 1, ^ + i-^ corresponding the position 
function of an equilateral triangle in the plane; we do this so each vertex is 
joined to precisely one of the same colour. We do not know which of the 
other roots lie in the spectrum since the computer program used for this 
example fails to solve the complete set of equations in a reasonable time. 
There remain two conjugate complex roots of the 7-polynomial. 

For the spectral values 7 = 1, we can compute the corresponding curva- 
ture as given by Definition 17.11 and Proposition 17.41 to obtain (in radians) 
5 = 27r at each vertex, to give a total curvature of 127r. 

As we calculated above, the bipartite graph ^^33 has 7-polynomial given 

by 

97^ - 357^ + 437 - 17 = (7 - 1)2(97 - 17) , 

In particular, these two graphs of constant degree three on six vertices cannot 
be isomorphic. 

9. An ELEMENTARY UNIVERSE 

Our objective is to construct an elementary universe populated entirely by 
graphs from which geometry and dynamics emerge. The universe is based on 
a binary relation between objects: are they connected by an edge or not? It 
is only this relation that matters; the nature of the objects being irrelevant. 
Our perspective is that, out of the graphs that are so formed, further implicit 
structure is present, given by the geometric spectrum and the corresponding 
fields. This implicit structure comes into play when graphs correlate. Thus 
we describe ways in which graphs can interact and so dynamics, that is 
change, appears. With a suitable definition of time, this change can be 
ordered to give a universe endowed with local geometry and time. 

The initial data for our universe is a graph T made up of a finite number 
of connected components Ti, . . . , Tk- In addition to the binary relation 
between vertices (whether or not they are connected by an edge), there is 
an additional relation between them: whether or not they belong to the 
same connected component of F. A priori there is no reason to give greater 
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emphasis to the property that two vertices be connected by an edge, rather 
than that they are not so connected. The representation by drawing an edge 
just gives a convenient way to visuahze the relation. 

A particle is a connected component of T. A state of the particle is an 
equivalence class of solutions (v?,7) to equation ([1]) on F^, where ip ^ Xtp + fi 
for A, /i G C and where we require that 7 < 1 at each vertex. A member of 
an equivalence class will be called a representative state and we shall write 
[if] for the equivalence class determined by the state if. A state for which 7 
is constant will be called an isostate. We allow point particles, consisting of 
a single vertex with all elements of C as representative states. 

A particle is not deemed to be in any state, but carries with it, its ensemble 
of states. By analogy with quantum mechanics, when two particles correlate, 
each falls into a particular state; that is, states are chosen with a certain 
probability. Isostates are favoured for empirical reasons relating to energy 
that we discuss at the end of this section. After correlation, further states 
may be present in the combined graph, permitting new correlations with 
other graphs that were not possible prior to correlation. An evolution of the 
universe is a sequence F 1— t- F' 1— t- F" 1— t- • • • of graphs, whereby a subsequent 
graph is obtained from the previous one by specific rules to be defined. The 
order of the sequence is dictated by the rules. There are three changes in 
our universe that we now specify: correlation between particles; internal 
mutation of a particle; separation of a particle into two or more particles. 

(i) Correlation. A correlation between two particles Fi = {Vi,Ei) and 
Fi = {V2,E2) is a new particle Fi ★ F2 = {Vi U V2,E), where we require 
EiU E2 C E, i.e. vertices are preserved and the edge set is increased. For a 
correlation to occur, we require Fi be in a state [(^1], F2 be in a state [(p2\ and 
F be in a state [ip] such that (p\y^ G [^pi] and ^p\v-^ G [^2]- After correlation, 
the new particle carries its ensemble of states and is not considered to be in 
any particular state. An example of correlation is shown in the figure below. 




-\/3i 2 - \/3i 
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Fig. 1. An example of correlation: the two particles at the top correlate to form a 

single particle in an isostate. 

We have drawn the resulting particle as an invariant framework in the 
plane, so the two central vertices appear superimposed at the point 1; but 
in fact they are distinct unconnected vertices. This correlation is to be 
favoured, since from two non-isostates, the outcome is an isostate. 

Separation. A separation of a particle is a dissociation of T into two par- 
ticles Fi = {Vi,Ei), T2 = {V2,E2) with F a correlation of Fi and F2. For 
separation to occur, we require F be in a state Fi be in a state [ifi] 
and F2 be in a state [(^92] with if\v^ G [ipi] and ip\v2 £ [^2]- An example of 
separation is given by the last example of Section [8l when the two concentric 
triangles connected by edges as indicated, falls into the state corresponding 
to the spectral value 7 = 1. The edges joining vertices on which the field 
has a common value are removed to give two disjoint triangles. 




Fig. 2. An example of separation 



Mutation. A mutation of a particle F is a change F = {V, i?) — )• S = (W, F), 
where S is a new particle with V = W. For mutation to occur, we require F 
be in a state [ip] , S be in a state [tp] with G [■0] . Collapsing is a particular 
example of mutation provided it does not disconnect the particle, whereby 
we remove edges that connect vertices on which if takes on the same value. 
If collapsing disconnects the particle, it falls into the category of separation. 
An example of mutation is illustrated in the following figure. 




Fig. 3. The particle on the left mutates into an isostate by the addition of an edge. 
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Another illustration is given by the particle of Figure 1, which, after 
mutation, produces the 1-skeleton of a cube. 



Fig. 4. Mutation occurs when two edges "flip" to connect the middle outer 
vertices to different central vertices. 

This change of state embeds the particle in Euclidean space as a 3- 
dimensional object. The particles before and after mutation both correspond 
to holomorphic states and so we view this mutation as neutral. 

Example 9.1. An evolution of a simple universe is as follows. First take 
the universe F consisting of two copies of the left-hand particle of Figure 
3. These then mutate to form F consisting of two copies of the right-hand 
particle of Figure 3. A correlation then occurs as in Figure 1 to form F". 
Finally a mutation takes place as in Figure 4 to form F'". We view this 
evolution as irreversible, in the sense that if we begin with the 1-skeleton 
of the cube, the probability that it fall into a state ip which has identical 
values on two diagonally opposite vertices to enable the reciprocal mutation 
of Figure 4, would be negligible. However, this is speculative, since we have 
not given a rule for deciding the probability of transition. 

A desirable objective would be to produce a complex universe from a sim- 
ple initial state. One way to accomplish this is to suppose the existence of 
virtual point particles that are susceptible to correlate with existing parti- 
cles. We now explore this possibility in more detail. 

Consider a graph F = (V, E) together with a function 99 : 1/ — )• C not 
necessarily a solution to ([1]). For ease of representation, suppose that each 
vertex x G ^ be placed at its corresponding position (^(x) in the complex 
plane. We fix our attention on a particular vertex, xq say, which by transla- 
tion, we suppose placed at the origin. Suppose xq has k neighbours placed 
at zi, . . . , Zfc. We now wish to add a new vertex placed at w and to join it 
to Xq in such a way as to satisfy ([T]) at xq. Specifically, we wish to consider 
the locus of points w which can be placed in this way. Since the mapping 



is in general holomorphic in w, we expect a 1-parameter family of values of 
w for which the right-hand side is real. 





(44) 



zi'^ + ■■■ + Zk^ + 



w I—)- 




A CLASS OF QUADRATIC DIFFERENCE EQUATIONS ON A FINITE GRAPH 67 



Example 9.2. Consider the graph on three vertices as indicated below. It 
may be that the extremal vertices placed at re'^ and 1 are joined to other 
vertices, but for the moment we are just interested in satisfying ([1]) at the 
origin. 




If we set w = u + iv, then it is a routine computation to show that the 
identity 7(1 + re'^ + w)"^ = 1 + r^e^'^ + w"^ has 7 real if and only if the 
following algebraic equation of degree three in u and v is satisfied: 

(nr sin — v{l + r cos 9)){u'^ + v'^) + rsin9{u'^ — v^) — 2uvrcos6 

(45) +r sin 9{1 — — 2r cos 9)u + (1 + r cos 6 + cos 6 + r"^ cos 26)v 

+r(l — r^) sin0 = . 

There are two cases when the solution set can be explicitly written down: 
(i) 9 = 7r/2. Equation now becomes: 

(46) {ru-v){u'^ +v'^)+r{u^ -v'^)+r{l-r'^)u+{l-r'^)v + r{l-r'^) = 0. 

For each r 7^ 0, 1, this is a smooth curve except at the singular point (n, v) = 
(— 1, — r). For the case r = 2, the curve is as indicated in the graph below. 
When r = 1 (so the original graph is holomorphic at the origin), this gives 
the algebraic set: 

{u - v){u^ + v'^ + u + v) = , 
consisting of the union of the line u = v and the circle (u+|)^ + (t>+^)^ = \ . 

A variant on the above procedure is to consider two particles F and S and 
to add a new vertex x which correlates with both particles, joining it to xq in 
F and yo in S, say. In order to correlate, suppose F falls into state [^p] and S 
falls into state [■(/']• Now we require that ([T]) be satisfied at the new vertex x. 
From Section[3l this is the case if and only if \(p{x) — ip{xo)\ = \'4>{y) — 4'{yo)\-, 
for representative states. This can be further generalized by creating a new 
vertex x and attempting to join several vertices to x. Such correlations can 
lead to discrete phenomena when we combine the various constraints. That 
is, the various loci determined by the real solutions to ()44p will in general 
intersect in a discrete set of points. 

The possibility that point particles may attach themselves to existing 
(more complex) particles, can lead to duplication and eventually a complex 
universe. The following sequence of correlation, mutation and separation 
gives an example of duplication. 
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We begin with a triangle on the left-hand side. Three isolated vertices then 

attach themselves in a symmetric way. This must be done to preserve the 

property that equation ([1]) remain satisfied. Symmetry is prefered since this 

leads to an isostate. In fact, if the vertices of the left-hand triangle have 

1 P\ 

representative field values 0, 1, 2 + 2 i) then the point particle connected to 
should have representative field value | + -^i, with the other point particles 
similarly assigned values to give an isostate with 7 = 1 (so the figure is 
misleading if we view the field as the position function, but avoids crossing 
edges). A mutation now occurs whereby the new vertices are joined by edges 
in the way shown. This produces an non-isostate particle with 7 = 7/9 at 
the new vertices (with 7 still equal to 1 at the original vertices), which then 
falls into an isostate given by the example on the left-hand side of Figure 
2, with spectral value 7 = 1, once more. Finally, separation occurs as in 
Figure 2. 

Further examples of how particles may correlate are given in Appendix 
1X1 This principally concerns unstable double cones which correlate to form 
a rich array of stable geometric structures. 
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We now wish to discuss how we may associate an energy to a state and the 
role of isostates. A natural quantity that occurs, which we may call energy, 
is given in the notation of Section [6] by = trace Z^Z (see equation 

m)- 

Recall that at each vertex, the vector Z picks out the configured star 
that projects to a state normalized to be zero at the vertex in question; 
furthermore, it has the minimum Frobenius norm amongst solutions to the 
system (|36|) . When = 3, it is precisely this solution that satisfies the 
constraint (j37p . In dimension > 3, the solution X = Z + Y which 
satisfies the constraint is no longer in general minimizing, whereas Z is still 
minimizing amongst solutions to ()36p . In dimension A^ = 2, at a vertex 
of degree 2 as discussed prior to Theorem 16.11 the solution to the lifting 
problem is completely determined. If we are to regard this energy as a 
criterion for stability in our universe, then this may provide reasons why 
dimension three should emerge as a favoured dimension. In defining energy, 
we also need to take into account the normalizing freedom. 

Definition 9.3. Let T = (V, E) he a connected graph endowed with a non- 
constant solution ip to ^ with ^{x) < 1 for all x G V (a state). Consider a 
particular vertex x £ V and let yi, . . . , y„ be the neighbours of x. Suppose ip 
is not constant on the star with internal vertex x. Set zi = ip{yi) — ^{x) for 
i = 1, . . . n, so that ()27p is satisfied for some 7. Then for N > 3, provided 
7 < 1, we define the energy of (p at x to be the quantity: 

£{^, X) := = |n - 7(3 - 27)^^ | ; 

22e\M 2(1-7) I iZelzer J 

for n = N = 2, we define the energy by 

£{<^, x) = 1 + cos 9 , 

where 6 is the exterior angle at x. If (p is constant on the star with internal 
vertex x, we take the energy to be zero. The total energy is defined to be the 
sum over the vertices of the energies at each vertex: £{(p) = Y2xeV ^i'^'-^^- 
There are various ways to normalize; we have chosen to divide by the 
average length of the complex numbers Z£. This is most convenient when 
n = N = 2, when the energy has the above concise expression. For A^ > 3, 
the formula in the definition is deduced from (|39p : 

I = trace Z^Z = p + a{l — — U2^) , 

and the expressions for ui and U2 given by (|30p . where we recall that a = 
7/>/(l — 7). For n = 2, we have ui^ + = 1, so that it is reasonable to 
replace \\Z\\^ above by p. From Sectional specifically equation ([U]), whenn = 
2, any solution to ([TJ with 7 < 1, necessarily corresponds to a configuration 
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in the plane in which the two edges connecting the neighbours of x have the 
same length, r say. But now if we refer to the figure above equation ([6]), then 

£{^,x) = '- = l- ^. ^. |l-e-^r = l + cosg. 

r 2(cos & — I) 

In this case, we can allow states for which 7 = 1 which correspond to 9 = zbvr; 
these are characterized as having zero energy. It is also important to note 
that the energy doesn't depend on the sign of the exterior angle. We now 
wish to look more closely at the case when n = = 2 in order to understand 
the importance of isostates. 

Consider a framework in the plane whose underlying graph is cyclic with 
K edges. Suppose the length of each edge is identical, so that the position 
function if of the framework defines a solution to ([1]). We are interested in 
critical configurations for the energy: 

K 

£{ip) = K + ^cosej , 
i=i 

where 9j is the exterior angle at vertex xj. In Appendix [Bl we show that 
the regular configurations, that is the isostates, are critical. 

Up to normalization, a three sided figure is completely determined and 
corresponds to an isostate with 7 = 2/3. Four sided figures are determined 
up to two branches by one of the exterior angles 6. The two branches 
occur depending upon whether we choose +9 or —9 for the opposite exterior 
angle. In the former case, the total energy is 4, whatever the exterior angle 
6 G (0, vr); in the latter case it is 2 + 2cos0. The two branches coalesce 
when 9 = 0. The absolute minimum £ = occurs when the four-sided 
figure is completely folded up, that is when all exterior angles are itvr. This 
corresponds to an isostate with 7=1. 

As K increases, the situation becomes more complicated. Up to normal- 
ization, the configuration space of the framework is parametrized by — 3 
exterior angles, however there are various branches that can occur which 
become more numerous as K gets larger. When K is even, the absolute 
minimum of £ is again zero and occurs when the figure folds up so all edges 
are superimposed and the exterior angles are all zbvr. When K is odd, the 
framework can no longer fold up in this way, but as we discuss in Appendix 
IHl the evidence suggests that the absolute minimum of £ occurs when the 
figure folds up as best it can, that is with all exterior angles as close to ±7r 
as possible, in a regular configuration. For example, in the case of a five 
sided figure, one can easily check that the regular pentagon gives a local 
maximum for £, whereas the regular star pentagon gives a local minimum. 
These are precisely the isostates for a five sided figure (see Example 13. 9p . 
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Time. Time is an ordering on a sequence of universes: (r,r ,...). The 
ordering must be compatible with the rules for change. Thus r^-'"^^) must 
derive from T^^^ by correlation, separation, mutation, or correlation with 
virtual point particles. There are two ways to decide such an ordering: 

(i) The rules for change: these may determine an irreversible process, 
such as that given in Example 19.11 The order (F, F', F", F'") is determined 
by the irreversibility of F" — t- F'" . 

(ii) A statistical parameter. The thermal time hypothesis has been devel- 
oped by Connes and Rovelli [12] . This is based on the Tomita flow associated 
to a von Neumann algebra. In quantum field theory, the appropriate von 
Neumann algebra is the closure of the algebra of observables. Then, given 
a state of a system over this algebra there is always a flow by which the 
state evolves and we may call this the "flow of time" (see also |26] §5.5.1). 
In our context, we don't have an obvious von Neumann algebra that we can 
exploit. However, there are various parameters that we may consider. 

Graph entropy is a well-know concept based on a probability distribution 
associated to the vertices pOj. Entropy is of course intimately related to 
the second law of thermodynamics. Intrinsic curvature is also a natural 
parameter that occurs in smooth Riemannian geometry and curvature flow 
provides a way by which a manifold may evolve into one of uniform structure. 

Curvature can provide a measure of local concentrations of structure. In 
general, our various notions of curvature described in Section [7] depend on 
the fleld ip satisfying ([T]). However, given a particle in a particular state, we 
could envisage processes whereby the graph could evolve to uniformize the 
curvature - say the vertex curvature, edge curvature or scalar curvature. 

A simple curvature which depends only on the combinatorial structure is 
given, for a graph F = (V, E) with degree function n : 1/ — )• N, by the func- 
tion 77,(2;) — 2. That this can be considered as a measure of curvature appears 
to have first been suggested in [32]. In [2], an algorithm is given whereby 
the quantity X]a;ey("'(^) ~ 2) 

may be minimized subject to ^^^ixav '^(■^) 
maining constant, by a process of sliding edges. This procedure preserves 
connectedness and may be viewed as a discrete analogue of the scalar cur- 
vature flow in Riemannian geometry. This leads to the parameter 

2\E\ 

as a possible measure of thermal time. In Example 19.11 the passage from 
F — 7- F increases t, whereas the passage from F — t- F decreases t, but we 
don't preclude local increases in t. Indeed, time should be a statistically 
dominant parameter that appears at a macroscopic level. 
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Appendix A. Invariant structures 

For our purposes, an elementary invariant framework is a framework cor- 
responding to the 1-skeleton of a either a regular polytope, or the invariant 
double cones of Section [SI possibly with the two apexes connected with an 
edge. An invariant structure is an invariant framework made up of ele- 
mentary components, the components connected by edges in an appropriate 
way. Motivated by the last section, we also require that the corresponding 
function 7 be constant. Thus, invariant structures are objects that may 
populate our elementary universe. 

Consider an invariant star in M.^ with internal vertex located at the origin 
and with n external vertices. Let x G M'^ be the centre of mass of the 
external vertices. Suppose that x ^ 0. Then we call the ray through the 
origin generated by x the axis of the star. Let 6 > denote the distance of 
the centre of mass from the origin along this axis. 

Lemma A.l. The addition of a new external vertex at any point other than 
—nb along the axis of the star produces a new invariant star. Furthemore, if 
7 denotes the invariant of the original star and x G M is the position along 
the axis of the new vertex, then the new star invariant is given by 

~ ^ {n + l){x^ + nb'^-f) 
^ ' ^ {x + nbY 

Proof. Without loss of generality, we may suppose that the centre of mass 
of the star lies along the i/Ar-axis. In particular, if -ui, . . . ,i;„ denote the 
external vertices, then 

n 

^ U£ = nbeN ■ 

i=\ 

We now add a new vertex at the point xe/v, for some x G M. Thus the new 
star matrix is given by 

{Vx\ ■ ■ ■ \Vn\xeK!) . 

As usual, let A = (ajk) be an arbitrary orthogonal transformation of and 
let P : — C be the projection P{yi, . . . , yN) = yi+iy2- Set Z£ = PoA{v£) 
for i = 1, . . . ,n and Zn+i = P o A(ae/v). Then 

N 

Zi = ^(aij + ia2j)v£j , Zn+i = x{aiN + ia2Ar) . 
i=i 

Furthermore, Yll=i ^t- — nb{aiN + ia2Af), so that 

2 

= jnb'^{am + ia2Nf , 
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where 7 is the invariant of the original star. We require that there is a real 
number 7 such that 



2 n 



n 



e=i / i=l 
But this is uniquely given by (HZ]). □ 

We note that as x approaches —nb, then 7 becomes arbitrary large. In- 
deed, when X = —nb, then we have harmonicity at the internal vertex of the 
new star, so that 7 is not well-defined in this case. 

Invariant structures now arise by connecting invariant frameworks with 
edges in an appropriate way. We can use the invariant double cones of 
Section [SJ as well as regular polytopes to produce new structures. There 
are various ways in which this can be done; as we don't have an exhaustive 
classification, we will consider some examples of geometric interest. 

Consider first the case of a double cone on a regular polygon. The corre- 
sponding framework in satisfies ([T]) with 7 in general having a different 
value at the apexes to the value at the lateral vertices. We extend the double 
cone by adding a new edge to each apex along the axis of the cone, attaching 
new double cones to each of these to produce an infinite family of double 
cones along a common axis. Suppose the double cone is aligned along the 
xa-axis. Let x denote a variable along this axis which is zero at the topmost 
apex and which increases towards the centre of mass of the vertex figure. 
Place a new vertex at position x along the axis. Then from ()47p and ()23p . 
the new value of 7 at the apex becomes: 

_{n + l){x'^ + n siiP ^ - f ) 

7apex — ; : 2^;^ ■ 

(X + n sm —j'^ 

We require this to equal the lateral value of 7 given by ()22p . This determines 
the quadratic equation in x: 

2x2 - n cos ^ 2 - 2 cos ^ + cos ^ 
^ ^ (x + nsin^)2" (2-cosf)2 

When n = 3, this has two distinct roots given by x = \/3 and x = — \/3/4. 
When we take the root x = \/3, then x lies precisely at the bottom vertex 
and the edge joining x to the top vertex connects the two apexes. In fact, 
since the triangle is a complete graph on three vertices, we have recovered the 
case of Corollarv l5.4l with n = 3. However, we can continue to add additional 
copies of the double cone to obtain a curious structure, as follows. 

We begin with one double cone Ci and attach it to a new one C2 so the 
top apex of Ci is joined to the bottom apex of C2 at a distance \/3 along 
the axis; thus the bottom apex of C2 is situated at the bottom apex of Ci. 
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The two cones are therefore superimposed, with a new edge running down 
the centre. The value of 7 at the top apex of Ci is equal to that at the 
bottom apex of C2, which is equal to 7iat- 

However, when we join the cones, we are at liberty to perform an arbitrary 
rotation of C2 with respect to Ci without affecting invariance. If we perform 
a relative rotation through an irrational multiple of 27r and perform the 
same relative rotation on joining a new cone C3 to C2 and so on, we obtain 
countably many double cones with common central axis. The closure of 
this set consists of the solid region enclosed by two (round) double cones, 
together with the segment A of the central axis from ^/S to 0. We note 
that even though A constitutes a common edge joining successive cones, as 
a framework, each edge is distinct. 

For n = 4, equation (j48p has no real solutions. This also turns out 
to be the case for n = 5, however for n > 6, there are two distinct real 
solutions, which enable us to connect double cones with edges along the axis 
of symmetry to obtain invariant structures. For example, for n = 6, we 
obtain the quadratic: 



Taking the positive root, on connecting successive cones we can perform a 
rotation to obtain interlaced frameworks. 

As a final construction, we reconsider double cones on regular polytopes 
and attach new double cones both laterally and vertically, interpolating the 
two distances in such a way that 7 is constant and the resulting structures 
are invariant. In order to do this, we need to be able to periodically position 
the regular polytopes in R^~^ in an appropriate way. For example, for 
regular polygons, this amounts to finding an appropriate periodic tiling in 
the plane. Three examples are illustrated below. 



118x2-48\/3 2;-27 = 



with roots, one positive, one negative, given by 



X = 



24^/3 ±^ 
118 




\ 
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In the first figure, we consider the square, bisect its vertex figure at each 
vertex to form a complementary octagon and so tile the plane with squares 
and octagons. For the second figure, we take a triangle, bisect its vertex 
figure at each vertex to form a complementary dodecagon (12-sided figure). 
Finally, the complementary figure of a hexagon is the hexagon itself. Our 
aim is to form an infinite structure in M^, by first taking the double cone 
on each of these polygons (square, triangle and hexagon); then taking an 
identical copy and attaching it to the first by vertical edges of appropriate 
length which join the apexes, so forming succesive layers. 

The first observation, is that the height of the double cone as given by 
Theorem 15.21 is not affected by the addition of a new edge and vertex bisect- 
ing the vertex figure. Indeed, using the notation of the proof of Theorem 
15.21 and taking the more general situation of that theorem of a double cone 
on a regular polytope, the star matrix is adjusted to 





( 


Vl 


V2 






d 







s = 




c 


C 




C 


a 


a 


i) 




\ 













b 


-b 





where d denotes the distance of the new vertex along the axis of symmetry 
(measured away from the centre of the vertex figure). In addition to the 
projections zi, . . . , Zn, Zn+i, Zn+2, we now have an additional vertex which 
projects to 

Zn+3 = -d{aiN + ia27v) • 
The affect of this is to give the new sums: 
n+3 

Z£ = {nc + 2o - d){aiN + ia2Af) , 

i=i 

and 

n+3 

ze^ = {nc^ + 2a^ + d^ - p){aiN + ia2Nf + (26^ - p){ai^N+i + ia2,7V+i)^ ■ 

Thus, as before, invariance requires b = ^ p/2 and the new value of 7iat is 
given by 

^ _ (n + 3)(nc2 + 2a2 + -p) 
^''^^ ~ (nc + 2a - df • 

If we take as an example the double cone on a regular triangle whose 
vertices are at a distance 1 from its centre, then we have 
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If we impose the tiling given by triangles and dodecagons, then d = ^/3, so 
that 

^ _ 40 

" (5 - ^3)2 ■ 

On the other hand, we can calculate 7 at the apex, after we add on another 
edge according to Lemma lA.ll To do this we apply ([26|) and ([17|) to give 

_ {m + l){x'^ + mb"^ — p-p) 8x^ + 6 

^''P'' " {x + mhY " (V2x- 37372)2 ' 

where we recall m = 3 is the cardinality of the regular polytope V (in this 
case the triangle) . Equating 7iat and 7apex yields the two (real) roots of the 
quadratic: 

4(9 - 5\/3)j;2 + 60\/3 X - 93 - 15\/3 . 

We now stack layers at either of these distances in order to obtain an invari- 
ant structure. 

One can generalize this procedure to higher dimension, whenever we can 
fill out the corresponding Euclidean space in an appropriate way. This can 
be done in any dimension for the hypercube and the cross-polytope, where 
we attach edges along the axes of symmetry of the various vertex figures. 
In M^, we can also take the tetrahedron with vertices placed at the points 
given by (fT6|l . Note that these lie at those vertices of a cube (of edge length 
2) diagonally opposite across each face. We then fill out with identical 
cubes, whose vertices are placed on a lattice with odd integer components. 
Begin by placing a tetrahedron in one of these cubes; then the axes of 
symmetry of the vertex figures cross adjacent cubes along diagonals which 
connect vertices which are opposite with respect to the centre of the cube. 
We obtain the 3-dimensional analogue of the tiling of the plane by squares 
and octagons illustrated above. By taking double cones on these tetrahedra 
and connecting them at the apexes with edges of appropriate length, we 
then obtain an invariant structure in M^. We omit the detailed calculations 
of the edge lengths, which proceed as in the example above. 



Appendix B. Some planar trigonometry 

In this section we parametrize the configuration space of the framework 
corresponding to an M-sided planar polygonal figure with edges of common 
length, and show that the regular figures are extremal with respect to the 
energy functional 8 defined in Section [9l 

Consider such a framework with M > 5 with edges of unit length labeled 
as in the figure below. The first bar has endpoints and 1 in the complex 
plane, then the next M — 3 exterior angles are labeled by cJi, (T2, . . . , (Tm-'a S 
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[— 7r,7r] and the last three exterior angles by 61,62,03. We define the quan- 
tities 

g = |]_ -|- gif^l ei{'^i+'^2) _|_ . . . _|_ gi(o-i+o-2H l-o-M-3)| 

J- — \1 -\- e''^^ + e^^'^'-^'^^^ + • • • + — i-o'm-4)| 

^ _ jgi""! -|- gKo'l+O's) _|_ . . . _|_ gi(o'l+(T2H |-CTAf_3)| 

= h -|- gi<^2 _|_ gi(o'2+<T3) _|_ . . . _j_ gi(o-2+o-2H ho-M-s)! 




Proposition B.l. The angles o"i, . . . , (Tm-3 parametrize a closed polygonal 
bar framework if and only if s < 2, in which case cos 62 is determined and 
provided s ^ 0, cos Oi and cos O3 have a two-fold ambiguity; the energy 

(49) 

£ := M + (^Yljt^^ cos a-j^ + cos 61 + cos 62 + cos 63 = 

+^^^ {E^I^[sin(ai + (72 + • • • + a-j) + sin(aj + aj+i + ■■■ + (Jm-s)]} 

where the choice of sign of the square root corresponds to two possible con- 
figurations of the framework. If s = 0, there are infinitely many possible 
configurations given by 62 = ±7r with 9i arbitrary. 

We postpone the proof of this proposition until the end of this section. 

Note that the sum of the exterior angles may jump from 27r to 47r with 
a continuous deformation of the framework, however, £ varies continuously, 
as indicated in the sketch below for M = 5. 




For the framework on the left, the sum of the exterior angles is 2tt, whereas 
for the right-hand figure, it is Att. At the point of transition, we have ai = 
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CJ2 = 27r/3. Then both s and the curly bracket (which equals 2(sin^ + 
sin ^)) vanish in the expression for £. However, the singularity is removable 
and £ = 3 — \/3 is well-defined and continuous at this point. 

Corollary B.2. The regular polygon and star polygons represent extrema 
for the functional E. 

Proof. The angles ai are not all on an equal footing, so in order to determine 
critical points, we establish a recursive formula on the derivatives. First note 
that 

s = |1 + e''^^ + e'*-'^^'*''^^^ + • • • + e'^'^^"'''^^'' — '"°"a/-3)| 

= i/{(l + coscTi + cos(cri + (T2) H h cos((Ti + (T2 H + iTm-s))^ 

+(sin cji + sin((Ji + (T2) H h sin(cTi + (T2 H h ctm-s))^} 

M-3 Af-3 

= - 2 + 2 ^ Yl + ^^^^+1 + • • • + ^i)} 

k=l j=k 

where the latter equality follows from using cos ai cos((Ti+cr2)+sin(Ti sin(o"i+ 
(72) = cos (T2 and so on. Therefore, for each £ = 1, . . . , M — 3, 

ds 1 ^ 

k=i j=e 

and we have the recursive formula 

-1 , M-3 



ds ds 



+ - ^ sin((Tfc H h fJ^-i) - - ^ sin((T£ H h 



5cj£ 9crf_i s f-' s 

k=l j=£ 

where, for i = 1, we set the first two terms on the right-hand side equal to 
zero. By the same reasoning, for i = 1, . . . , M — 4, 

or dr 1 ^ . , N 1 ■ / n 

— = K - > , sm((Tfc H h cr<?_i) > sm((T£ H h ctj) , 

fe=l j=£ 

with dr/daM-3 = 0, and for ^ = 2, ... M - 3, 

5f 5t 1 4-4 , . 1 , . 

— = H 7 > , sm((Tfc H h (T£_i) - 7 > ^ sm((T£ H h cTj) , 

k=2 j=t 

with dt/dai = 0. From (jl9|) . we now obtain for ^ = 2, . . . , M — 3, 
(50) 

— = h smcT^^i - sm(T£ - - sm(cri H h fi^^i) 

acr£ oa£^i 2 



+1 sm{ae H h cTM-3) - ^ ^s''" [cos{ai H h <t^-i) - cos(fT£ H h <tm-3 
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where we have set 

M-3 

T := ^ [sm(cJi H h o-j) + sm{aj H h ctm-s)] ■ 

i=i 

We now claim that under the hypothesis that ah angles = 9 are equal 
with M9 = 2kTr, for some integer k satisfying < k < [M/2], then all 
derivatives dS/dai vanish. 

The case 6 = zbvr can only occur when M is even, in which case it rep- 
resents an absolute minimum for £ and is certainly critical; so henceforth, 
suppose that 9 ^ ztyr. With all angles equal, we must take the positive sign 
for the square root in (09]); indeed, this is necessary to obtain M — M cos9 
for the value of £. The following identities are useful: 
1 + gie + e2ie ^ . . . ^ g(M-i)ie ^ g 

sm{j9) = - sin(M - 2)9 - sm{M - 1)0 = sin 0(1 + 2 cos 9) 
similarly ^^l^^ cos{j9) = — cos 0(1 + 2 cos 9) 



T = 2 sin 0(1 + 2 cos ( 



and 



2 cos — r = t = |l + 2 cos ( 
2 ' 



On noting that dt/dai = 0, for the derivative with respect to we obtain 

d£ . r dr A ds . , . „, 

— = -sm(Ti + -- --^==— > sm(j0 



A/ -3 



I cos(M - 3)0 + ^ cos(j0)} 



2s 



0, 



where the last equality follows after substitution and routine calculations. 
Now apply recurrence on dS/dai, so suppose that d£/da£^i = 0. Then 
from ()50p we have, 



d£ I 1 COS0 

— = — sin(^- 1)0- -sin(2 + £)0 + [sin ^0 + sin(£ + 1)01 

oae 2 2 1 + cos 

Sin 

But now multiplication through by 2(1 + cos 0) shows that this vanishes. By 
recurrence, all derivatives d£/da£ vanish for £ = 1, . . . , M — 3. □ 

Included in the list of "star polygons" are ones that correspond to several 
circuits of a regular polygon. For example, when M = 6 and in the notation 
of the above proof = 2tt/3, the solution corresponds to two circuits of a 
regular triangle as illustrated in Example 13.81 Such cases occur when M has 
non-trivial factors. When M is even, we always have the solution = zt-zr. 
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where the polygonal chain "folds up" to cover the segment from to 1. Then 
the energy is zero and is at an absolute minimum. We conjecture that the 
regular convex polygon given hy 9 = 2tt/M is a global maximum for £ and 
for M odd, the star polygon given hj 9 = 2[M /2]'k /M is a global minimum; 
but we don't have proofs of this in general. Also, we do not know if there 
are other critical configurations for £ other than the regular ones. 

When M = 5, there are just two possible regular configurations. For 
<7i = o"2 = 27r/5, we have £ = 5(3 + \/5)/4 and for ui = (T2 = 47r/5, 
we have £ = 5(3 - \/5)/4. It is now possible to calculate explicitly the 
second derivatives of £ at these critical points and to verify directly that 
these configurations correspond respectively to a local maximum and local 
minimum of £. However, this approach for general M presents formidable 
computational difficulties. 

Our introduction of the functional £ suggests the intriguing possibility 
of setting up a gradient fiow on the configuration space of positions of the 
framework parametrized by the cr^, whereby the figure evolves into one of the 
regular configurations (depending upon whether we take the gradient flow 
or its inverse). The combinatorial solution to such problems of regularizing 
polygonal frameworks, has been derived in the references [J^ and [9]. But 
we are not aware of a functional analytic approach. An evolution from the 
pentagon to the star pentagon with reflectional symmetry, is obtained by 
setting o"! = (72 and allowing ai to vary continuously from 27r/5 to 47r/5. 
Two of the positions of this evolution are illustrated in the figure above. 



Proof of Proposition. Without loss of generality we can take the length of 

2 



each edge to be unity. Since s = 2 cos cos 02 is completely determined by 



the data ai, . . . , (Jm-z and is given by 

cos U2 = 



2 

even though 02 is only determined up to sign. We now find expressions 
for 01 and ^3. Suppose s 7^ 0. Note that this implies that r ^ as well. 
Let a be the angle indicated in the figure, being the interior angle of the 
triangle with adjacent edges of lengths r, s and opposite edge of length 1. 
Set A = 9i + Then application of the sine rule gives the equalities: 

sin(j4 — a) sin A 
sin a = = . 

s r 

One can eliminate the terms involving a from this expression to obtain 

— — 1 



cos^ 



2s 

Once more this is unambiguously defined, even though A itself has a sign 
ambiguity. This is underlined by the expression for cot A, which has the 
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ambiguity attached: 

cot^ 



— — 1 



^^4^2 _ (^2 _ ^2 _ 1)2 

However, 

cos 9^ = cos(j4 ) = (1 + cot A cot — ) sin A sin — 

^2 2 2 

from which we deduce that 

1 + cot ^ cot % 

cos 6i = 



;i+COtM)(l+COt2^) 
= -^{^4 - sV4s2 - (r2 - s2 _ 1)2 + 5(^2 _ ^2 _ ^ 

But now there is a reflectional symmetry of notation by traversing the poly- 
gon in the opposite sense, so we may interchange aj with ^ with t 
and 9i with ^3, to obtain 

cos 03 = ^{\/4 - sV4s2 - (t^ - - 1)^ + s{t'^ -s^- 1)} ■ 

Simphfication of the various terms now gives the formula of the proposition. 
When s = 0, then 62 = ±7r and we can pivot the two superimposed edges 
that result about a common vertex, giving infinitely many solutions obtained 
by varying 9i arbitrarily. □ 

Appendix C. The linearized and weak forms of the equations 

We develop the relevant functional analytic framework on a finite graph 
in order to deduce weak forms of the equations ([T|). This provides heuristic 
arguments as to why a pair (F, (p) consisting of a finite connected graph and 
a solution (/? to ([1]) should model an elementary particle. 

Let r = {V, E) be a finite graph. For x £ V define the tangent space to T 
at X to be the set of oriented edges with base point x : T^T = {xy : y ~ x}. 
Define the tangent bundle to T to be the union: TV = Ux^yTx^- Then 
a 1-form on F is a map oj : TV — ?• C such that u;{xy) = —uj{yx). To a 
function : 1/ — t- C, we can naturally associate a 1-form, the derivate dip, 
by difixlj) = ip{y) - ip{x). 

For two 1-forms to, r], define their pointwise symmetric product at x £ V 

by 

{uj,Tj)x = ^uj{xy)ri{xy) , 
and their (global) symmetric product by 

('^> ^) = X] '^(^)^(^) = ^ XI X] ^(.xy)r]{xy) . 
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Note that in the first sum the 1-forms act on unoriented edges so that only 
their product is well-defined; the factor of one half occurs in the second sum, 
since there, unoriented edges are counted twice. 

For functions (p,ip : V ^ C, define their (global) symmetric product by 

(9?, V') = ^ n{x)ip{x)ij{x) , 

where n{x) is the degree of vertex x. 

The above definitions are the complex symmetric analogues of standard 
products that arise in functional analytic theory on a graph; in the latter 

situation they are replaced by Hermitian products rather than symmetric 

products 

Given a function ^ : y — )• C and a 1-form uj : TV — )• C, we can define a 
new 1-form by 

{(,oj){xy) = ^{C{x) + ^{y))uj{xy) . 

Then it is easily checked that 

d{ipip) = ifdip + ipdip . 

Given a 1-form w, define its co-deriviative d*(jj to be the function which 
at each vertex x G ^ is given by 

Then for a function 99 : y — )■ C, 

d*d(/3 = ]— V dipixy) = ^ - = -A99 . 

n(x) ^-^ n\x) ^-^ 

^ ' y^x ^ ' y^x 

Lemma C.l. Let Lfjtp : V ^ C be functions and uj : TT — >• C a 1-form. 

Then the following formulae hold: 

(i) {dip,uj) = {ip,d*^) ; 

(ii) (A(p,V) = -(d(p,dV); 

(iii) (Av9,^) = (v?,AV'); 

(iv) d*{ipu}){x) = ip{x)d*uj{x) - ^:^{d<p,uj)x for each x eV. 
Proof. To prove (i), we notice that for x ~ y, the sum 

xeV y^x 

contributes {ip{y) — (p{x))oj{xy) = —Lp{x)uj{xy) — {p{y)uj{y~x) (the term being 
symmetric in x and y), which equates to the corresponding terms in the sum 

{(f, d*uj) = - ^ ip{x) u}{xy) . 
xeV y^x 
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The identity (ii) now follows from the fact that d*d<p = —A(p. Identity (iii) 
follows from (ii), by symmetry. Finally 

d*{ipuj){x) = ^ V((^a;)(x1/) 

nix) ^-^ 



ip{x)d*u){x) - ^{^{y) - <fix))uj{xy) 



= ip{x)d*u{x) - ^^^^{dip,uj)x , 
which gives (iv). □ 

With the above notation and formulae established, we can give a weak 
form of equation ([1]) . 

Proposition C.2. The equation dTJ holds if and only if 
(51) (A(^,^7A(^) - 2(d(/9,<d(/?) = 0, 

for any function ^ : y — ?• C. 

Proof. Let ^ : ^ — t- C. Then if ([1]) holds, we have: 

^ n(x)e(x) {j{x){Aip{x)f - difixf) = , 
xev 

equivalently 

{A^,C^Aip)-{d^^O = 0, 
where we recall that dip'^{x) = Y^y^^^i^iy) ~ — (d^^, d(/?)a;. But we 

claim that (d(^^,^) = 2{dip,S,dip). Indeed, 

(dv.^o = E ^(^)^(^) ( - ^(^))' ) ' 

xeV \y^x / 

which for each x ~ y, contributes the term {n{x)^{x) + n{y)^{y)){ip{y) — 
ip{x))'^. But precisely this term occurs with in 

2(dv9, n^d^) = IY1 E(^(^)^(^) + n{y)^{y)){ip{y) - ip{x)f . 
x£V y^x 

(on noting the symmetry of the expression to be summed on the right-hand 
side). 

Conversely, if we fix a vertex x £ V and consider the function ^ : 1/ — )■ C 
given by ^{x) = 1 and ^{y) = for all y ^ x, then ()5ip gives 'y{x)Aip{x)'^ — 
dip'^ix) = 0, so that holds. □ 

On taking the function ^ to be identically equal to 1, we obtain the 
following consequence. 
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Corollary C.3. Let ip : V ^ C be a solution to equation ([TJ with 7 constant. 
Then 

hA(A(p) + 2nA(p + -{dn,d(f) , (f) = 0. 
n 

This corollary suggests heuristic arguments as to why we might consider 
a pair (F, (p) consisting of a connected graph endowed with a solution ip to 
([TJ with 7 constant as a particle with mass inversely proportional to I7I; 
in the case when 7 = 0, we will view the pair as representing a massless 
particle. 

Firstly, we do not admit any fixed background with respect to define 
parameters of equations: the particle in a sense creates its own background, 
so we view an equation of the form {T'{ip),ip) as appropriate, where V is 
some (discrete) differential operator. In the case when n is constant, we now 
note the relation between the operator V{(p) = ^A(Aip) + 2nAip and the 
operator on the left-hand side of the time-independent Schrodinger equation 
on a fixed smooth background: 

(-^V^ + V{x))^ = E^{x), 
\ 2m J 

when ip is identified with A(p. The case of mass zero (when 7 = 0) is justified 
in some detail in [3]. 

In order to define the linearized equation, we consider a family {(ft} of 
functions such that ipQ = ip solves ([1]) with 7 independent of t. On writing 
^(x) = ^^^f^\t=o, we obtain the equation linear in ^: 

(52) 7(x)Av9(x)AC(x) = {dipix),dC{x)), . 

In the case when 7 = 0, that is ip is holomorphic, this takes the particularly 
simple form: 

(53) (d(^(x),de(x)),. = 0, 

for all x €zV. In all case, we see that ^ = Xip+fi solves the linearized equation 
(A,/x G C constant); this reflects the normalisation freedom ip 1— t- Xip + fi. 

The natural class of mappings between graphs which preserve equation ([1]) 
are the so-called holomorphic mappings. These were introduced for simple 
graphs under the name semi-conformal mappings by Urakawa \32\ I33j. as 
the class of maps which preserve local harmonic functions (harmonic at a 
vertex). The notion was later extended to non-simple graphs by Baker and 
Norine [11[5]) who used the term holomorphic mapping. In [3j, it was shown 
that the holomorphic mappings are precisely the class of mappings which 
preserve local holomorphic functions, as we have defined them by ([2]). 

The definition requires that we restrict to mappings of the vertices that re- 
spect the relation of adjacency. Thus we define a mapping f : T = {V, E) ^ 



A CLASS OF QUADRATIC DIFFERENCE EQUATIONS ON A FINITE GRAPH 85 

S = (W, F) between graphs as a mapping of the vertices such that x ^ y 
impUes either f{x) = f{y) or f{x) ~ /(?/)• 

Definition C.4. Let f : T = {V,E) — t- S = (W,F) be a mapping between 
graphs. Then f is holomorphic if there exists a function X : V ^ N such 
that for all X £ V and for all z' ^ z = f{x), we have 

\{x) = A(x, z) = 'i^{x G y : x' ~ X, f{x) = z'} , 

is independent of the choice of z' ; we set X{x) = if f{x') = z for all x' ~ x. 
Call A the dilation of /. 

Proposition C.5. Let f : T = (y,E) ^ = {W,F) be a holomorphic 
mapping between graphs of dilation A : ^ — >• N. Suppose il) : W ^ C 
satisfies the equation 

MAV)' = (d^)2, 

for some /u : — )• M. Then for each x & V such that X{x) / 0, the function 
(p = ijj o f satisfies ^ at x with 

^ n{x)fi{f{x)) ^ 
A(x)m(/(x)) ' 

where m{f{x)) denotes the degree of the vertex f{x) G W . 

Proof. Let f :T = {V, E) ^ T, = (W, F) be a holomorphic mapping between 
graphs of dilation A : y — N. Let x £ V and set z = f{x). Then 



Ii{z) 



m{z) 

Since / is holomorphic 

^ [(^ o f){x') - o f){x)] = A(x) i^Piz') - i:{z)) . 

Suppose that A(x) 7^ 0. Then 

^[(V^o/)(x')-(^o/)(x)f = X{x)Y^{^{^)-^{z)f 

Y^{^{:^)-i,{z))\ 



m{z) 
X{x)m{z) 



from which the formula follows. If on the other hand A(x) = 0, then f{x') = 
f{x) for all x' ~ X and both sides of ([1]) vanish. □ 



86 paul baird 

References 

[1] W. W. Adams and P. Loustaunau, An Introduction to Grobner Bases. Graduate 

Studies in Math. Vol. 3, Amer. Math. Soc, 1996. 
[2] P. Baird and M. Tiba, An algorithm to prescribe the configuration of a finite graph, 

arXi^ 1002:1192 vl. 

[3] P. Baird and M. Wehbe, Twistor theory on a finite graph, Comm. Math. Phys., 304, 

No. 2 (2011), 499-511. 
[4] M. Baker and S. Norine, Riemann-Roch and Abel-Jacobi theory on a finite graph, 

Advances in Math. 215, No. 2 (2007), 766-788. 
[5] M. Baker and S. Norine, Harmonic morphisms and hyperelliptic curves. Int. Math. 

Res. Not. (2009), 2914-2955. 
[6] S. Barre, Real and discrete holomorphy: introduction to an algebraic approach, J. 

Math. Pures Appl. 87 (2007), 495-513. 
[7] B. Buchberger, Ein Algorithmus zum Auffinden der Basiselemente des Resklassen- 

ringes nach einem nulldimensionalen Polynomideal, Ph. D. Thesis, Inst. University 

of Innsbruck, Innsbruck, Austria, 1965. 
[8] F. R. K. Chung, Spectral graph theory, CBMS Regional Conference Series in Mathe- 
matics, vol. 92, Washington DC, 1997. 
[9] R. Connelly, E. D. Demaine and G. Rote, Straightening polygonal arcs and con- 

vexifying polygonal cycles. Discrete and Computational Geometry 30, No. 2 (2003), 

205-239. 

[10] R. Connelly, T. Jordan and W. J. Whitely, Generic global rigidity of body-bar frame- 
works, EGRES. TR-2009-13, www.cs.elte.hu/egres/. 

[11] R. Connelly and W. J. Whitely, Global rigidity: the effect of coning. Discrete and 
Computational Geometry 43, No. 4 (2010), 717-735. 

[12] A. Connes and C. Rovelli, Von Neumann algebra automorphisms and time- 
thermodynamics relation in general covariant quantum theories. Classical and Quan- 
tum Gravity 11 (1994), 2899-2917. 

[13] H. S. M. Coxeter, Regular Polytopes, Third Edition. Dover Publ., New York, 1973. 

[14] R. Diestel, Graph Theory, Graduate Texts in Math., 3rd Ed., Springer- Verlag, Berlin, 
Heidelberg, 2005. 

[15] M. G. Eastwood and R. Penrose, Drawing with complex numbers, Math. Intelligencer 
22 (2000), 8-13. 

[16] R. Ehrensborg, The Perles-Shephard identity for non-convex polytopes, 2007: 

www.scientificcommons.org/42923842 
[17] C. F. Gauss, Werke, Zweiter Band, Koniglichen Gesellschaft der Wissenschaften, 

Gottingen 1876. 

[18] M. Gromov, Metric Structures for Riemannian and Non-Riemannian Spaces, 2nd 

Edition. Birkhauser, Boston, Basel, Berlin, 2007. 
[19] B. Griinbaum and G. C. Shephard, Descartes' theorem in n dimensions, Enseign. 

Math (2) 37 (1991), 11-15. 
[20] J. Korner, Bounds and information theory, SIAM Journal on Algorithms and Discrete 

Math., 7 (1986), 560-570. 
[21] J. Milnor, The Schldfli differential equality, John Milnor Collected Papers, Vol. 1, 

Geometry. Publish of Perish Inc., Houston Texas, 1994, pp. 281-295. 
[22] J. Murakami, On the volume formulas for a spherical tetrahedron, larXiv:1011 .2584vl. 

Nov. 2010. 



A CLASS OF QUADRATIC DIFFERENCE EQUATIONS ON A FINITE GRAPH 87 



[23] R. Penrose, A generalized inverse for matrices, Proc. Cambs. Phil. Soc. 51 (1955), 
406-413. 

[24] R. Penrose, Angular momentum: an approach to combinatorial space-time, in Quan- 
tum Theory and Beyond, ed. T. Bastin, Cambridge Univ. Press 1971, pp. 151-180. 

[25] C. Rovelli, Relational quantum mechanics, International Journal of Theoretical 
Physics 35 (1996), 1637-1678. 

[26] C. Rovelli, Quantum Gravity. Cambridge Monographs on Mathematical Physics, 
Cambridge University Press, 2004. 

[27] G. T. Sallee, Stretching chords of space curves, Geometriae Dedicata 2 (1973), 311- 
315. 

[28] L. Schlafli, On the multiple integral f" dx dy ■ ■ ■ dz, whose limits are pi — aix + biy + 

h /n« > 0, p2 > 0, . . . ,p„ > 0, and + y^ -\ z^ <1, Quart. J. Math 2 (1858), 

269-301. 

[29] G. C. Shcphard, Angle deficiences of convex polytopes, J. London Math. Soc. 43 

(1968), 325-336. 

[30] I. Streinu, A combinatorial approach to planar non-colliding robot arm motion plan- 
ning, Proc. 41st IEEE Annual Symposium on Foundations of Computer Science 
(FOGS), 443453, 2000. 

[31] B. Sturmfels, What is a Grobner basis? Notices of the AMS, 52 No. 11 (2005), 2-3. 

[32] H. Urakawa, A discrete analogue of the harmonic morphism, in Harmonic morphisms, 
harmonic maps and related topics. Research Notes in Mathematics, vol 413 (ed. C. 
Anand, P. Baird, E. Loubeau and J. C. Wood), 97-108. Chapman and Hall/CRC 
Boca Raton, 2000. 

[33] H. Urakawa, A discrete analogue of the harmonic morphism and Green kernel com- 
parison theorems, Glasgow Math. J. 42 (2000), 319-334. 

Departement de Mathematiques, Universite de Bretagne Occidentale, 6 AV. 
Victor Le Gorgeu - CS 93837, 29238 Brest Cedex, France 
E-mail address: Paul.Baird@imiv-brest.fr 



